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1. Introduction 

In all fields of physics homogeneous systems have the simplest properties and thus 
play a particular role. This also holds with respect to phase transitions and critical 
phenomena. The general behaviour, the universal features near second-order transitions 
and the universality classes in this case are well-known [1]. On the other hand, all real 
systems are inhomogeneous in one way or another. This may affect their properties in 
different ways. Enhanced couplings in a finite region of a ferromagnetic system will result 
in a locally stronger magnetic order. For the case of competing interactions the effects 
of an inhomogeneity may be more complicated in detail but will still be of a local and 
quantitative nature. If, however, the perturbation has sufficient extent, the character of 
the phase transition and the critical behaviour may change. 

One example for this situation is well-known: a system with a free surface. To obtain 
it from a homogeneous system one has to cut an infinite number of bonds. The case of 
planar surfaces was first treated for the two-dimensional Ising model [2, 3] and subsequently 
studied in great detail [4, 5]. It was found that, connected with the surface, there is a set 
of critical exponents with values different from those in the bulk. The modified critical 
behaviour is seen in a boundary layer which has a width of the order of the correlation 
length. 

The planar surface, however, is only the simplest example showing such an effect. 
A number of other inhomogeneous systems display similar features. It is then of interest 
to collect and present these cases, to see relations between them and to discuss the basic 
aspects and mechanisms. This is the aim of the following review. 

The systems to be discussed are mainly classical spin systems with short-range 
ferromagnetic interactions and, occasionally, walks and polymers. The inhomogeneities 
have a regular, non-random nature and the examples fall into two groups. One consists of 
systems in various geometrical forms with free boundaries described by algebraic curves, 
the other contains models with defect lines or extended inhomogeneities with algebraically 
varying strength. In both cases unusual and nonuniversal critical behaviour is found if the 
perturbation is sufficiently effective. Then the local exponents depend continuously on a 
parameter or the behaviour changes altogether from power laws to streched exponentials. 
In the second group of systems also local ordering at or above the bulk transition 
temperature may occur. The necessary condition for these effects always involves scaling 
dimensions of the undisturbed system. 

There are homogeneous systems which also show such features. An example is the 
eight- vertex model [6]. In this case, however, all exponents vary in a similar way and 
one can obtain universal values by measuring the temperature via the correlation length. 
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The inhomogeneous systems considered here have a richer structure since local and bulk 
quantities appear simultaneously. In those cases where the effects are related to the 
geometry they also are of a rather general nature since non-universal behaviour occurs 
for any simple isotropic system and in all dimensions. 

The mentioned phenomena can often already be seen in a mean field treatment 
which is therefore given in various places. More generally, renormalization and scaling 
considerations allow to classify the perturbations and the way they act. There also is 
a rather general method to calculate the change of critical exponents via perturbation 
theory. A majority of the examples are two-dimensional systems and for them two other 
methods are available. Conformal invariance makes detailed predictions at the critical 
point and has been used in various cases although its validity in inhomogeneous systems 
is not necessarily guaranteed. Exact calculations for various Ising systems, however, have 
always confirmed conformal results. Moreover they allow to obtain a complete picture 
by giving results for all temperatures. The main tool here are transfer matrices and the 
two-dimensional problems are thereby related to certain quantum spin chains. All these 
topics are treated in the various sections and in three appendices. The latter contain some 
background material and some more technical aspects. 

Certain types of inhomogeneity will not be considered here although they may lead to 
interesting effects. Random systems, by their very nature, demand special methods and 
are still under investigation with respect to their critical properties. For reviews we refer to 
[7, 8] . In the case of wetting phenomena one is concerned with the properties of interfaces 
which are created by appropriate boundary conditions and which may be influenced by 
inhomogeneities. This is also a field in itself and has already been reviewed [9-11]. Finally, 
layered systems with a periodic structure [12-14] will also be omitted. 

Some of the effects to be discussed should be observable in experiments but the article 
deals with the theoretical aspects. Nevertheless, it is not aimed at the specialist. The 
notation has thus been chosen as simple and coherent as possible and therefore differs 
from the usual one in some places. 
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2. Planar surfaces 



At a free surface missing bonds reduce the order and, as mentioned before, the critical 
behaviour is modified locally. The loss of full translational invariance leads to a decay 
of correlations which is different parallel and perpendicular to the surface. In dimensions 
higher than two, enhanced surface couplings may induce surface order above the bulk 
critical temperature. The local critical behaviour can also be influenced by the surface 
geometry or through the introduction of long-range surface induced perturbations. These 
last two points will be considered in the next sections. 

The easiest approach to the problem at hand makes use of mean field theory which 
often brings a good qualitative description of the phenomena and becomes exact above 
the upper critical dimension where the fluctuation effects are negligible. Since a detailed 
account of planar surface critical behaviour within this frame can be found elsewhere [15, 4] 
we limit ourselves to a brief study of the bulk transition and the ordinary surface transition 
for illustrative purpose. This will serve as an introduction to the mean field treatment of 
the other, less usual, inhomogeneous systems which are considered in the next sections. 

There also exists a number of recent reviews dealing with surface critical behaviour 
beyond mean field theory [4, 5, 16, 17]. As a consequence, only selected results concerning 
the structure of correlations, critical profiles and surface exponents are discussed in this 
section. 



In Landau mean field theory the system is treated in a continuum description. The 
total free energy F of a sytem with a volume (V) limited by a surface (S) is written as the 
sum of bulk and surface contributions which are functionals of the order parameter m(r) 
which is nonvanishing in the ordered phase 



Let us consider for simplicity an Ising system with a scalar order parameter m and a 
free energy which in zero external field is even in m, i.e. symmetric under order parameter 
reversal. Near a second order transition, the order parameter is small and the bulk free 
energy density fb[m] is written as an expansion in the order parameter and its gradient, 
limited to the following terms 



2.1. Mean field theory 




(2.1) 




(2.2) 
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The second one, with C > 0, gives the extra energy associated with a spatial variation of 
the order parameter, A ~ T — T c measures the deviation from the critical temperature T c , 
B is a positive constant to ensure stability below T c and h is the bulk external field. 
In the same way the surface free energy density is written phenomenologically as 

(2.3) 



f 777 ^ 

f 8 [ m }=f s [0] + -C—-h 8 m 



where m is the value of the order parameter on (S). The quantity A, with the dimension of 
a length, is called the extrapolation length and can be deduced from the microscopic surface 
and bulk interactions through a mean field treatment of the microscopic Hamiltonian of 
the system [4]. 

The equilibrium value of the order parameter m(r) minimizes the free energy in (2.1). 
It may be obtained through a variational calculation by looking for the first order change 
of the free energy, 8F[m], associated with a deviation 5m(r) of the order parameter from 
its equilibrium value. Using (2.1-3), one obtains 



5F[m] = / [CVm -V5m + (Am + Bm 3 - h)6m] dV 

J(V) 



+ 



(S) 



—m - h t 



(2.4) 



5m dS. 



The first term in the volume integral may be rewritten as 

CV/ii • VStit =V-(C<5mVm) - CV 2 m 5m 



(2.5) 



and the contribution to (2.4) of the first term on the right can be transformed into a surface 
integral through Gauss theorem. Then 



5F[m] = / [-CV 2 m + Am + Bm 3 - h] 5m dV 

J(V) 



<(S) 



_ m h s 
-n ■Vm+- — — 

A C 



(2.6) 



5m dS 



where n is a unit vector normal to the surface and pointing inside the system. 

At equilibrium the first order variation of the free energy vanishes. The volume part 
gives the Ginzburg-Landau equation 



CV 2 m(r) = Am{r) + Bm 3 {r) - h 



(2.7) 



governing the bulk equilibrium behaviour whereas the surface part provides the boundary 
condition 



n • Vm(r) 



m(r) h 



s 

C 



(2.8) 
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In the bulk the l.h.s. in (2.7) vanishes and the zero-field magnetization is given by 

_({-A/B)W~t 1 /* T<T C 
mb ~\ T>T C (2 - 9) 

where t ~| A \ is the reduced temperature. The spontaneous magnetization vanishes at T c 
as a power of the reduced temperature with a mean field bulk exponent (3=1/2. 
The connected part of the order parameter correlation function 

SfTli V ) 

< m{r)m{r') > c =< m{r)m(r') > — < m(r) >< m(r') >= (2.10) 

Oft/if* J 

can be obtained through the introduction of a varying bulk external field h(r) in (2.2). 
This simply amounts to replacing h in (2.7) by h(r). Taking a derivative with respect to 
h(r') and using the definition (2.10) leads to 

[-CV* + A + 3Bm 2 (r)] < m{r)m(r') > c = 8(r - r'). (2.11) 

For constant m(r), which is the case in a homogeneous system below T c or generally 
above T c , the correlation function is therefore proportional to the Green function which 
satisfies 

(-V 2 r + ^jG(r,r') = 8(r-r'). (2.12) 



Here 



T > T c 

£=< v^- (2.13) 

is the bulk correlation length diverging at the critical point with a mean field exponent 
v = 1/2 which is the same in both phases. Hence the critical G is the Green function of 
the Laplace equation which establishes a link to electrostatic problems. 

For a flat surface at y = without external surface field, the boundary condition (2.8) 
gives 

dy A 

When the surface interactions are not larger than in the bulk, due to the missing couplings 
at the boundary, the magnetization always increases from the surface into the bulk and 
the extrapolation length is positive. The surface transition is then driven by the bulk and 
there is an ordinary surface transition at T c . Since A <^ £ near the bulk critical point, 
the magnetization profile extrapolates to zero near the surface and one may use Dirichlet 
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boundary conditions, i.e. m = on the surface, instead of (2.14). With rh = m/nib and 
h = 0, (2.7) and (2.9) lead to 

d 2 m A A S 2 o A. o „ , . . 

= —m + —m h m = — — (m —m). (2-15) 

Multiplying both sides by 2dm/dy and integrating, one obtains 

dm 1 — m 2 

*7 = — <2 ' 16) 

where the integration constant has been chosen to give a vanishing slope at infinity and 
the form of the correlation length in (2.13) has been used. The solution of (2.16) with 
rh = 1 at infinity then is 

y 

m(y) = m b tanh— . (2.17) 

In the vicinity of the surface, m ~ t y, so that the magnetization grows linearly with the 
distance to the surface. Its temperature dependence is different from the bulk one and 
involves a new exponent (3 S = 1 which is the mean field surface exponent at the ordinary 
transition. 

The form of the correlation function for the ordinary transition is given by the solution 
of (2.12) which satisfies Dirichlet boundary conditions at the surface. To obtain this 
function for a half-space, one can start from a (i-dimensional slab limited by two surfaces 
at y = and y = L. With r — r' = r\\ + {y — y')n and taking advantage of translational 
invariance in the d—1 transverse directions, one may rewrite G(r\\,y,y') as a Fourier 
expansion with components G^{y,y') satisfying 

(-^2+^G k {y,y') =8{y-y>) (2.18) 

where k = \/k 2 + £ -2 . The solution can be written as an eigenfunction expansion 



2 sm(mry/L) sm(mry'/L) 



G k {y,y) = Jj ^ k2 - n2n2/L2 

n=l ' 



(2.19) 



For a semi-infinite system L — > oo and the sum over n can be replaced by an integral 
which is evaluated using the method of residues so that finally [15] 

G(r h y,y>) = J d d ^k ^ [e^'l - e"^] . (2.20) 

The two parts of G are each correlation functions of the infinite system, one between r 
and r' and the other between r and the image point of r' relative to the boundary. At 
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the critical point they decay with power (d — 2). From this the surface-bulk correlations 
follow by taking r\\ = 0, y' fixed and y ^> y' , leading to the asymptotic decay G(y) ~ 
y-(x+x s ) ^ ^-(d-i) Here x and x s are the bulk and surface (ordinary) scaling dimensions 
of the magnetization (Appendix A.l). The bulk behaviour is obtained taking y and y' 1 
with | y — y' |>> 1. Then the decay exponent is 2x = d — 2 so that 

x = — x s = — (ordinary transition). (2.21) 

The scaling relations (A4, A7) (3 = vx = l/2 and f3 s = i>x s = l are satisfied by the mean field 
exponents at the upper critical dimension which, in this case, is d c = 4. 

Finally, consider the critical profile of the magnetization when its value at the surface 
is fixed at m(0) = mo- It follows from (2.7) with A = and h = so that 

Following the same steps as above, the integration gives 

1/2 



m{y) = m 
with the asymptotic behaviour 



(2.23) 



m{y)~y-\ (2.24) 

The magnetization decays as a power of the distance to the surface, a behaviour which is 
also obtained beyond mean field theory as discussed in the next section. 

If the surface couplings are sufficiently enhanced and/or a surface field is present, A 
becomes negative and the surface orders at a temperature which is higher than the bulk 
critical one. Then bulk ordering at T c induces some singularity in the surface behaviour. 
The associated critical point is the extraordinary transition. Increasing A decreases the 
temperature of the surface transition until the surface and the extraordinary transition 
meet at the special transition. A detailed treatment can be found in references [4, 15]. 



2.2. Correlation functions and critical profiles 
Conformal methods (Appendix A. 2) can be used to make quite general statements 
about critical systems. They generalize, at the critical point, the covariance under global 
scale transformations which is at the basis of the renormalization group, by introducing 
local scale transformations with a varying dilatation factor b(r). Such local dilatations 
are realized via conformal transformations. In any dimension the conformal group is 
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constructed by including the inversion r' = r/r 2 besides the usual uniform transformations: 
translation, rotation and dilatation. 

Particularly useful in the study of surface properties is the special conformal 
transformation [18,16] 

PS = £ + a ( 2 - 25 ) 

which combines an inversion followed by a translation and a new inversion. A semi-infinite 
system with a flat surface containing the origin is invariant under such a transformation 
when the translation a is parallel to the surface. Using an infinitesimal translation, 
covariance under this transformation determines the form of the critical correlation 
functions [18] and allows a determination of boundary induced profiles. Although the 
transformation works in any dimension, only the two-dimensional situation is discussed 
below. 

Consider a critical system on a half-plane with a surface at y = 0. In Cartesian 
coordinates (2.25) translates into 

x' = x + e(x 2 -y 2 ) y' = y + 2exy (2.26) 

where the infinitesimal translation is a = (— e, 0). Using complex notations, z = x + iy, 
(2.25) is rewritten as 

z' = z + ez 2 (2.27) 



and the local dilatation factor is 

dz 



b(x) 



dz' 



1 - 2ex. (2.28) 



Let ip(r) be some local operator (energy density, order parameter) with bulk scaling 
dimension x, i.e. transforming as ip(r/b) = b x ip(r) under a global change of scale. Its 
two-point critical correlation function G(x\,X2,yi,y2) =< i>(xi,yi)i>(x2,y2) >, which is 
a function of x\ and x-i through u = x\ — xi due to translational invariance parallel to the 
surface, is transformed into 

G(u\ y[, y' 2 ) = b( Xl ) x b(x2) x G(u, yi ,y 2 ) (2.29) 

according to (A14) under (2.26). A first-order expansion in e leads to the differential 
equation 

dG dG dG 

{x\ -x\- y\ + yt) -JT- + 2x iyi — + 2x 2 y 2 — + 2x( Xl + x 2 ) G = 0. (2.30) 
du dyi ay 2 

10 



Inhomogeneous systems 

Since G depends on x±, x 2 only through u, all the terms involving a factor x\ + X2 sum 
up to zero and ones obtains two partial differential equations 



dG dG dG 

u— + yi- h y 2 ^ h 2xG = 

ou dy 1 dy 2 

The first one expresses the homogeneity of the correlation function 



(2.31) 



b 1 b 1 b 

or, with b = u and Q = yi/u, 



G[T^^)=b 2x G(u, yi ,y 2 ) (2.32) 



G(u, yi ,y 2 ) = «-^G(l,Ci,C2) = (^CiC 2 )- x S(Ci,C2). (2.33) 

Using this form, the second equation in (2.31) fixes the way spatial coordinates combine 
into a single scaling variable p in the scaling function S(p) so that, finally, 



- x 2 , yi, 2/2) = (2/12/2) x S 



2/12/2 



(xi - x 2 ) 2 + (2/1 - 1/2) 5 



(2.34) 



The behaviour of the scaling function for small or large values of the argument can be 
deduced from scaling considerations. For instance, with y\ = y 2 fixed and | x\ — x 2 | = | 
u |— > 00, the surface-surface correlations decay like | u \ ~ 2Xs where x s is the surface scaling 
dimension of ip. As a consequence 

H(p) ~ p Xs 0. (2.35) 

In two dimensions the scaling function itself satisfies a certain differential equation 
following from conformal invariance and has been determined explicitly for the Ising, Potts 
and O(N) models [18, 16, 19, 20]. 

The same method applies to the determination of critical profiles when tfj(r) is some 
quantity with nonvanishing average at the bulk critical point of a semi-infinite system. 
This may be either the energy density with free or fixed boundary conditions or the order 
parameter with fixed boundary conditions. Since in two dimensions the surface is one- 
dimensional, an enhancement of surface interactions is not sufficient to maintain surface 
order at the bulk critical point. Therefore the boundary variables must be fixed in order 
to have a non-trivial order parameter profile. 
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Due to translational invariance parallel to the surface, the critical profile depends only 
on y. Following the same steps as above, a differential equation for < ip(y) > is obtained, 
leading to the universal profile 

< ifj(y) >~ y~ x . (2.36) 

That (2.36) follows from conformal invariance was already noticed in a different way in 
[21, 22]. This algebraic form was originally deduced from simple scaling considerations 
[23]. In any number of dimensions, < ip(r±/b) >= b x < ip(r±) > where r± is the distance 
to the surface, so that with b = r± , 

< if>(r±) >~ rl x . (2.37) 

When ip is the order parameter, x = (3/u which is equal to 1 in mean field theory in 
agreement with Equation (2.24). 



2.3. Surface exponents 

Some heuristic arguments have been used [24] to obtain the values of surface exponents 
in arbitrary dimensions. They apply to the ordinary and extraordinary transitions for 
the energy density and to the extraordinary transition for the order parameter. In two 
dimensions the results are supported by independent conformal arguments. 

Consider a system inside a cube with L d interacting spins and free boundaries. 
Increasing the size of the system by SL through the addition of a surface layer, the partition 
function becomes 

Z L+5L = Trexp [-(3(H + 6H)} = Z L < exp(-[35H) > L (2.38) 

where 5H is the energy change associated with the extra layer. To leading order in L the 
free energy 

F L = L d f b + 0(L d - 1 ) (2.39) 

only involves the bulk free energy density Its variation under an infinitesimal change 
of the size is given by 

SF L = -/T 1 In < exp(-/35H) > L ~< 5H > L ~ dL^ 1 f b SL + 0(L d ~ 2 ) (2.40) 

so that 
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and the surface energy density e s on the l.h.s. scales near the bulk critical point like the 
bulk free energy density 

e s ~ ~ f b ~ t 2 " Q . (2.42) 

Thus its scaling dimension x s at the ordinary transition is equal to the dimensionality of 
the system according to the Josephson hyperscaling relation (A2). The same argument 
applies to the magnetization and energy densities with fixed boundary conditions, i.e. at 
the extraordinary transition. Then 

m s ~ e s ~ t 2 ~ a (2.43) 

so that both scaling dimensions are equal to the dimension d of the system. For the surface 
energy operator at the ordinary transition, this result was first obtained by Dietrich and 
Diehl in the O(N) model using renormalization group methods and the short-distance 
expansion [25]. 

In two dimensions, any analytic function w(z) on the complex plane provides a 
conformal transformation with a local dilatation factor b(z) = | dz/dw | (see Appendix A. 2). 
Such a transformation can be used to deduce surface exponents starting from the half- 
space critical profile. Under the conformal mapping w = L/n arcosh z the half-plane y>0 
is transformed into the half-strip (u > 0, < v < L, w = u + iv). The dilatation factor is 
b(z) = ir/LyJ sinh 2 (7ra/L) + sin 2 \irv / 'L) and y = sinh(7rtt/L) sm(iiv/L). The boundary- 
induced critical profile < ip(y) > (2.36) which has a bulk scaling dimension x, then 
transforms into [22] 

< ip(u, v) > ~ b(z) x y~ x 

in the half-strip. The profile in the new geometry can also be written as an expansion 
in terms of the eigenstates of the transfer operator T = exp(— H) where H is the strip 
Hamiltonian with appropriate boundary conditions (Appendix B.l). Then 

< ip(u, v)>=^2M n <n \ ij){v) | > exp [-(E n - E )u] (2.45) 

n 

where E is the ground-state energy of H and M n selects the eigenstates | n > compatible 
with the boundary conditions at u = 0. The smallest gap corresponding to a nonvanishing 
matrix element gives the surface scaling dimension x s of ip as in (A23). Since the large-w 
expansion of the r.h.s. of (2.44) involves only powers of exp (— 2ttu/L) the gaps in (2.45) 



+ sin 



2 



\ L 
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are multiples of 2rc/L. Assuming that M\ is non-nanishing, one may identify the surface 
scaling dimension as 

x s = 2 (2.46) 

in agreement with previous results in the special case d=2. 

The considerations presented above suppose the existence of a non- vanishing profile at 
the critical point. As a consequence they cannot be applied to the case of the magnetization 
at the ordinary surface transition. As mentioned in the preceding section, other conformal 
techniques can be used, which completely determine the critical correlation functions in 
the semi-infinite geometry. From these, the ordinary surface exponents given in Table 2.1 
for the g-state Potts model and the O(N) model have been identified [18, 16]. 



g-state Potts model (d = 2) 


Q 

x s 


12 3 

I I 2 

u 3 2 3 


4 
1 


O(N) model (d = 2) 


N 

x s 


-2-1 1 

i 13 5 1 
1 16 8 2 


2 
l 

4 



Table 2.1. Order parameter scaling dimension at the ordinary surface transition in the g-state 
Potts model and the O(N) vector model. 



3. Wedges, corners and cones 

Systems in the form of a wedge show new features and their study gives new insight into 
the influence of geometry on critical behaviour. They were first investigated by Cardy [26] 
using mean field theory in various dimensions and an e-expansion near d = 4. Subsequently 
the corner geometry in d = 2 was studied by accurate calculations, mainly on Ising models, 
and by conformal mapping. The main result is that the edge and corner exponents are 
functions of the opening angle 9 of the wedge. The planar surface and its exponents can 
be viewed as a simple special case. 

This non-universal behaviour can be related to the particular geometrical properties of 
a wedge, which does not contain a length parameter and thus is invariant under rescaling. 



14 



Inhomogeneous systems 



The opening angle is therefore a marginal variable in a renormalization transformation 
and may enter into the expressions for the exponents. The same will happen for all other 
scale-invariant figures like arbitrary pyramids or cones. Furthermore, the system does not 
have to fill the figure, it may just cover the surface. For such a situation there are even 
exact analytical results. 

3.1. Mean Geld theory 
Within a continuum approximation, in the mean-field approach one has to solve the 
same equations as in Section 2.1, but for a different geometry. The order parameter 
correlation function is proportional to the Green function (2.12) with Dirichlet boundary 
conditions for the ordinary transition. This function also appears in electrostatics [27] and 
in diffusion problems [28] . In the present context it was discussed by Cardy [26] for wedges 
in d dimensions. 

For a wedge, one proceeds as for a slab and performs a Fourier transformation in the 
d—2 directions parallel to the edge. In the remaining directions one uses polar coordinates 
p, a where a is the azimuthal angle in the wedge (0 < a < 9). The Fourier components 
Gfc then satisfy 

with k defined below (2.18). The solution can be written as an eigenfunction expansion 
and takes the form 

_ , . 2 >^ [°° dp a , , / a . fnixa\ . ( mva'\ /n _ 

G k (p,a;p,a)= q2^J q k 2 + ^2 ^^P) JAPP ) sin j sin J (3.2) 

where the J u are Bessel functions and v = mr/9. The continuous index p is related to 
the infinite extent of the system in the p-direction. Using this expression one can find, 
for example, the critical correlations parallel to the edge. Their asymptotic behaviour is 
determined by the contributions of n = 1 and small p. One obtains 

^ (3-3) 

with the edge exponent 

Xe{0) = ^ + |. (3.4) 

The decay thus is always faster than in a homogeneous system, where it would be given 
by the first term in (3.4) as found in (2.21). This is an effect of the boundary conditions. 
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In an electrostatic picture, G gives the potential at r due to a point charge at r' and the 
effect comes from the induced charge of opposite sign on the (metallic) boundaries. For 
9 = n/n, where n is an integer, it can be simulated by (2n — 1) image charges. As the 
wegde becomes narrower, the effect becomes stronger. 

To obtain the order parameter profile in the wedge one would have to solve the non- 
linear Ginzburg-Landau equation (2.7) with h = 0. For the critical behaviour, however, 
this is not necessary. The radial dependence has the general scaling form nibfip/Q an d 
since m is small, the spatial dependence follows from the linearized equation. This gives, 
for p < f 

m(p,a)~m 6 (|J sin (J^j (3.5) 
Inserting the temperature dependence of £ and m& then leads to the exponent 

This result is valid in all dimensions. At the upper critical dimension d c = 4 the scaling 
law j3 e = vx e as in (A7) is satisfied although j3 e and x e are non-universal. For 9 = tt the 
surface exponent j3 s = 1 is recovered. For smaller angles, the order parameter near the 
edge vanishes with zero slope at the critical temperature. This reflects the difficulty to 
maintain the ordered state in such a geometry. 

One can also treat a semi-infinite wedge where the third planar boundary is 
perpendicular to the edge. This system has a three-dimensional corner with two right 
angles and the third one equal to 9. The correlation function then follows from the image 
method as 

G(r,r') = G 00 {r,r')-G 00 {r,r") (3.7) 

where G^ is the result for the infinite wedge and r" is the image point of r' with respect 
to the third plane. For r' near the corner, the correlations parallel to the edge then decay 
with a power x c + x e = 1 + 2x e , from which the corner exponent x c (9) = 1 + x e (9) follows. 
If 9 = 7r/2, one obtains a cubic corner with three right angles and x c = 7/2. The exponent 
of the order parameter is calculated as above 

(3 C (9) = 1 + ^. (3.8) 

For small 9 this is the same law as in (3.6) but in general (3 C > (3 e so that the order near 
the corner vanishes faster with temperature than near the edge. In particular for the cubic 
corner one has j3 c = 2, i.e. a quadratic behaviour. 
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The general three-dimensional corner, for which the exponents depend on three angles, 
has not been treated. However, results can be given for a cylindrical cone of opening angle 
9 with respect to the axis. The Green function of this problem has been determined in 
other contexts [28, 29]. For r' on the axis and r > r' it has the form 

G(r,r') =J2 A » (y)^ Pm(cos«) (3.9) 

where again < a < 9. The Legendre function P M (cosa) has to vanish on the boundary 
(a = 9) which determines the allowed values fx = \i m , m = 1, 2, • • •. The asymptotic decay 
is again a power law with the exponent given by the smallest Hi = Hi(9). For the order 
parameter near the apex of the cone one obtains, in the same way as for the wedge 

a (0) = i[l + /ii(0)]. (3.10) 

The function Hi(9) is shown in reference [30]. For 9 —> it varies as 1/9 and the system 
has the same features as a narrow wedge. For 9 = n/2 one has \i\ = 1 and the result 
Pa = Pa = 1 for the planar surface is recovered. For 9 — > n the quantity ji\ vanishes 
logarithmically. Thus one obtains the bulk result (3 a = (3 = 1/2, although one is dealing 
with a system from which an half-infinite needle is cut out. 



3.2. Ising models 

A number of results exist for corners in two-dimensional Ising lattices with certain 
discrete values of the angle. The corner magnetization m c has been calculated for triangular 
and square lattices in three different ways : from star-triangle recursion equations [31], from 
the corner-corner correlation function [31-33] and from Baxter's corner transfer matrix [34]. 

In the recursion method, one starts from a finite lattice, e.g. in the shape of a triangle 
with fixed spins along one edge, and reduces it successively to smaller sizes. In this way 
the corner magnetization can be evaluated exactly, but the procedure has to be done 
numerically. In the second method one considers the correlation between spins at adjacent 
corners of a lattice in the form of a square. Asymptotically, it factors into the product of 
the two corner magnetizations for which an expression 

m c ~ 3.11) 

<0|S> v 1 

is found. Here the operator af refers to the corner spin, the states |0 > and |1 > are 
the eigenstates of the row transfer matrix (Appendix B.l) with largest and next-largest 
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eigenvalue, and \B > describes the state of free spins in the upper and lower boundary row. 
The quantity m c can then determined from the solution of a matrix or integral equation. 
The corner transfer matrix, finally, is already the partition function for a whole angular 
segment (see Appendix B.2). Thus it is, from a geometrical point of view, the ideal tool 
for treating the corner geometry. However, due to the edges of free spins one has to take 
a matrix element < B\T\B > and the calculation becomes relatively difficult. One has to 
solve a matrix equation also in this case. 
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Figure 3.1. Corner magnetization vs. temperature for isotropic Ising lattices with various 
opening angles. Except for the 60° case all results refer to square lattices, (a) edges along the 
bonds (b) edges along the diagonals. 



In Figure 3.1 results are shown for isotropic systems with various opening angles at 
the corner. The basic feature is that m c decreases with decreasing 9 for all non-zero 
temperatures. There is no general explicit result for m c (t,9) because even the transfer 
matrix calculations involve some numerics in the end. Based on the numerical results 
there is, however, a conjecture for a square lattice with a 90° corner and edges along the 
bonds [33] 

m c = 1 - -(cothif! - l)(cothK 2 - 1) (3.12) 

where Ki, K2 are the couplings in the two directions. The orientation of the edges and 
the precise location of the chosen spin both affect m c [32] but not the critical exponent for 
which the formula 

ft = £ (3.13) 
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was inferred from the data for isotropic systems with discrete angles [31]. Since v = 1 here, 
this is also the result for x c . Thus one finds a simple 1/9- dependence on the angle as in 
the mean field treatment. If the system is anisotropic equation (3.13) is not immediately 
valid. However, by a rescaling as in Appendix B.l, one can make the system effectively 
isotropic. The opening angle 9 thereby changes to an effective angle 6* e // 5 as illustrated in 
Figure B2, which need not be a simple fraction of ir. With the new angle, equation (3.13) 
is again satisfied. It was shown later that this ^-dependence is a consequence of conformal 
invariance (see Section 3.4). 

Ising models covering the surface of a pyramid have also been treated [31]. They 
are obtained by putting together k lattice segments, each with opening angle a, so that 
9 = ka 2-7T. The calculation in this case is much simpler since the system is closed upon 
itself and has no edges of free spins. One can then use corner transfer matrices and derive 
the apex magnetization m a from a simple trace formula, see (B18-19). The result takes 
the form 

9 



(21 - 1) e(t) 2 r 



(3.14) 



m a (t, 9) = J^tanh 
i=i 

where e(t) ~ l/ln(l/t) for t — > 0. Curves for m a (t) look qualitatively similar to those 
shown in Figure 3.1. The critical behaviour can be found by converting the product in 
(3.14) into the exponential of an integral. The logarithmic dependence of e on t then leads 
to a power-law behaviour of m a . From the way the angle 9 enters it follows that the apex 
exponent (3 a = x a is related to the bulk value via 

Pa = (3.15) 

This result also follows from conformal invariance. For a finite system at the bulk critical 
point, the variable t~ l in e, which is proportional to the correlation length, is replaced by 
the size L [35] so that m a varies as L~ Xa as in (A8). 



3.3. Other systems 

The corner geometry has also been investigated for the (9(A^)-model in the limit — > 0. 
In this case the spin correlation function is related to self-avoiding walks (SAW's) on the 
underlying lattice via [36] 

oo 

G(r, r') = RNw N(r, r') (3.16) 

N=0 
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where K is the coupling constant of the vector model and Wjv(r, r') denotes the number 
of SAW's with N steps, going from r to r'. By counting walks up to N ~ 25 for several 
angles, Guttmann and Torrie [37] found in two dimensions 

*c = |. (3.17) 

For wedges in three dimensions they found the edge index 

Xe = 0.5 + 0.847^. (3.18) 


So far this is the only result in d = 3 beyond mean field theory. 

One can also go the other way and derive properties of Wn from those of G, invoking 
conformal invariance [38, 39]. This was also done for polymers in corners [40]. For the walks 
one obtains, for example, the spatial moments of Wn- These were studied for systems in 
d = 2, 3 with excluded half-infinite lines, which are special cases of corners and cones, 
respectively [41, 39]. 



3.4. Conformal results 

In two dimensions the critical behaviour of a wedge can be related to that of a half- 
plane if the system shows conformal invariance. This is done through the conformal 
transformation w = z e ^ [18, 31]. 

The critical correlation function G(w, w') in the wedge is found by transforming (2.34) 
according to (A14). For w' near the corner and w far in the bulk, it varies as 

G ^ w ')~y^- ( 3 - 19 ) 

where x c is given by 

x c = ^x s . (3.20) 

Thus the corner index is related to the surface index via the geometrical factor rr/9 which 
also appears in the mapping. The results (3.13), (3.17) are special cases corresponding 
to the surface exponents x s = 1/2 and x s = 5/8 in the two models. This confirms that 
conformal invariance actually holds in these systems. By relating other exponents to x c 
and repeating this for different scaling operators one can then derive the whole set of wedge 
indices. 
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Figure 3.2. Two-dimensional corner geometry with two types of boundary conditions. 

In the same way a conical system is obtained by mapping a full plane onto a wedge with 
periodic boundary conditions via w = z e ^ 2ir . The correlation function in this geometry is 
then obtained from the one in the full plane via (A14). This leads to 

2n 

x a = —x (3.21) 

relating apex and bulk indices, as already found in the treatment of the Ising model, 
cf. (3.15). 

The mapping from the half-plane can also be used to determine the universal critical 
profiles, in the wedge geometry, which arise by fixing the variables along certain boundaries. 
Two situations are shown in Figure 3.2. For identical edges (Figure 3.2a) one finds [22] 

< (j)(p,a) >~ [psin^— JJ 0<a<9. (3.22) 

If one of the edges is free, the surface index x s also enters. In the half-plane this corresponds 
to a situation where one has free variables along half of the real axis and fixed ones along 
the other half. The profiles for such a case have so far only been determined for Ising 
and Potts models [42-44]. Taking these results one can also find the profiles for the case 
shown in Figure 3.2b. The necessary mapping is w = ( e / n with £ = a{yfz — l)/(y / i + 1). 
This transforms the upper z-plane onto a half-disc in the £-plane and onto a wedge with 
circular boundary of radius R in the w-plane. Near the corner, the result for the profile is 

X c — X TV j 

<(P(p,6/2)>~ p ° rxc /0 « J R. (3.23) 

For the analogous conical system the profile only depends on p and follows from the result 
for a disc [22]. The size dependence is then given by R~ Xa . 
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With complex mapping one can also study rounded corners. In agreement with general 
predictions [5] one finds that the rounding does not affect the corner exponents or the 
asymptotic form of G in (3.19). 

Conformal invariance also predicts universal contributions to the critical free energy 
which are connected with a corner or an apex [45] . For a finite two-dimensional system of 
characteristic size L, the free energy has the form 



where the terms written are the bulk and surface contributions, respectively. For systems 
with Euler number x = like cylinders or tori, the next term is a universal constant, 
related to the Casimir effect [46, 47]. For x 7^ 0, however, as is the case for simply 
connected domains, the next term is a logarithm in L. Each corner leads to a contribution 



where c is the conformal anomaly characterizing the universality class of the system. For 
the Gaussian model (c = 1) there is a close connection [48] between this term and the 
contribution of a corner to the eigenvalue spectrum of the Laplace operator as studied by 
Kac and others [49-51]. Thus the same geometrical factor appears in both cases. An apex 
of a cone gives an analogous contribution as (3.25), with the substitution tv/6 — > 2tt/9 in 
the bracket. 

Adding the contributions (3.25) for a system in the form of a rectangle gives AF = 
— (c/4) InL. For the Gaussian model with fixed boundary variables this can be checked by 
a direct calculation [52, 53]. This result has been used to discuss the shape dependence of 
the critical free energy [54]. For a polygon with a large number of edges, AF approaches 
the result 



valid for smooth boundary curves. Privman has given an interpretation of the logarithm 
in terms of an interplay of singular and non-singular contributions in the free energy [55]. 
According to this argument such terms can also arise from corners in three dimensions. 



F = f b L 2 + f s L + ■ ■ ■ 



(3.24) 




(3.25) 



AF = ——InL 

6 



(3.26) 
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4. Parabolic shapes 

A qualitative change in the critical behaviour occurs for systems which are 
asymptotically narrower than wedges or cones. This was observed in a study of two- 
dimensional systems bounded by algebraic curves of the form v = ±Cu a [56]. A typical 
case is the simple parabola (a = 1/2) and, generally, such shapes will be called parabolic. 
The wedge corresponds to the special case a = 1. Three typical situations are shown in 
Figure 4.1. 




(a) (b) 



Figure 4.1. Three types of two-dimensional systems bounded by simple parabolic curves. 

Using conformal and transfer matrix methods the following features were found. For 
a < 1 the critical behaviour is no longer characterized by the usual power laws but 
by stretched exponentials. Thus thermodyamic quantities show essential singularities. 
Furthermore there is non-universality since the amount of stretching is determined by the 
value of a. For a > 1, on the other hand, one finds normal surface critical behaviour. The 
difference between the two cases can be attributed to the way the dimensional parameter 
C behaves under renormalization for a < 1 and a > 1. The wedge geometry, a = 1, 
is thereby seen to be the borderline case. A corresponding distinction holds in higher 
dimensions. For intrinsically anisotropic systems an analogous classification holds but the 
marginal case is then found for a ^ 1. 

4.1. Mean Geld theory 
The simplest example is a two-dimensional system with parabolic boundary, v 2 = 
2pu + p 2 , as in Figure 4.1a. Its critical correlation function can be obtained by direct 
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calculation from (2.12) or using the transformation z = i cosh n^w/2p which maps the 
upper z-half-plane on the interior of a parabola in the w-plane. For points on the axis 
with u' fixed and u large, one finds 



G(u, u') ~ exp 



— 7T 



U 

2p 



1/2 



(4.1) 



The decay of G is therefore faster than in a wedge, but slower than in a strip where G 
varies as exp(— iru/L), if the width is L. The power w 1 / 2 in (4.1) is clearly related to the 
form of the boundary. The length scale is set by the parabola parameter p or, equivalently, 
by the quantity C = \/2p. 

In three dimensions, a similar result is obtained for a paraboloid of revolution [57]. 
Working in parabolic coordinates, one can express G in terms of Bessel functions [58]. On 
the axis it has the form, for u> u' 



G(u,u')=J2 A m(u')K 



m=l 



u 
2p 



1/2- 



(4.2) 



where \i m is the m-th zero of the Bessel function J . Therefore asymptotically 



G(u,u') 



u 



1/4 



exp 



u 



1/2- 



(4.3) 



with hi ~ 2.41. Here the power in front of the exponential is also related to the shape 
since it does not appear for a cylinder. As in Section 3.1, the result may be interpreted in 
electrostatic terms. The induced charge in the present case is more concentrated near the 
source point u' so that the decay of G is faster than in the conical geometry. 




Figure 4.2. Geometry of parabolic system near the tip. The shaded portion of size D along 
the axis is the region where the shape governs the critical behaviour. 
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There exists no exact calculation of the order parameter m{r) using (2.7). However, 
near the critical point one may integrate the linearized equation from the tip up to 
u = D ~ £, 2 /p, where the width of the system becomes comparable to £ and m(r) reaches 
its bulk value (Figure 4.2). In this region m increases exponentially in (u/2p) 1 ^ 2 and, for 
fixed u 

m ~ exp(— a£/p) (4-4) 
where a is a constant. The order therefore vanishes exponentially in t -1 / 2 at T c . 



4.2. Conformal and scaling results 
In two dimensions one can go beyond mean field theory and obtain results for arbitrary 
models which are conformally invariant. Complex mapping from the upper ^-half-plane 
allows to obtain the three types of shapes shown in Figure 4.1. For a< 1 (Figure 4.1a) it 
is given by z = icosh(w/2p) a . The critical correlation function is then found as in the case 
of a wedge. For points on the axis and large values of u [56] 



G(u,u') = A(u')-^exp 



nx. 



2C(1 -a) 



u 



l-a 



(4.5) 



The functional form of G thus varies continuously with the parameter a. The fall-off is 
fastest for a = which corresponds to a semi-infinite strip of width 2C. With increasing a 
the exponential becomes more and more stretched and approaches a power law for a — > 1, 
the wedge limit. For a = 1/2, the mean field (Gaussian) result (4.1) is reobtained by 
inserting the bulk and surface exponents x = 0, x s = 1 (2.21). One also notes that the 
three-dimensional result (4.3) has the same structure as (4.5) and the power 1/4 in (4.3) 
can be interpreted as ax, since x = | is the mean field exponent in three dimensions. 

By mapping the full z-plane one obtains systems with identical upper and lower 
boundaries ("paraboloids"). As in the passage from the wedge to the cone in Section 
3.4, one then has to replace ttx s by 2nx in (4.5). For the Ising and Potts models one may 
also find, as in Section 3.4, the complete order parameter profile at the critical point if 
the variables are fixed at the right end of the system. The size dependence of the order 
parameter near the tip is then given, for both geometries, by the same exponential which 
appears in (7, with u replaced by the size L. 

For a > 1 (Figure 4.1b) the mapping has to be changed into z = i(w s —p s ) l / s where 
s = (a — l)/a. Then on the axis, asymptotically, 



G («»« , )~^5 7 (4-6) 
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which is the result of the half-plane. Finally, if a parabolic piece with a < 1 is cut out of 
an infinite system (Figure 4.1c) one obtains 

G(uy)~-^r (4-7) 

with the corner exponent x c (2n) = x s /2 for a cut (see (3.16)). 

These results can be understood if one considers the behaviour of the boundary curves 
under a rescaling r — > r' = r /b. Their functional form, being a power law, is invariant but 
the parameter C changes according to 

C' = b a - x C. (4.8) 

For a > 1, C grows under renormalization and the boundary approaches a straight line. 
For a = 1, C is invariant and thus a marginal variable (the angle of Section 3). For 
a < 1, C decreases and the system approaches either a cut geometry or becomes locally 
one-dimensional. In the latter case there is a destruction of long-range order which leads 
to the particular features found above. The same considerations hold for paraboloids in 
three dimensions. A parabolic cylinder, on the other hand, renormalizes to a half-plane in 
the latter case which has long-range order. This order, however, appears only at a lower 
temperature and stretched exponentials at the original transition are still possible. 

The variable 1/C may be considered as a scaling field which plays the same role as the 
variable 1/L in finite-size scaling (Appendix A.l). It has dimension 1 — a and vanishes at 
the fixed point of the half-plane geometry. It appears in scaling relations like the one for 
the magnetization on the axis (cf (A31)) 

m (t, u, 1) = b-*m p *^p) (4.9) 

leading to the functional form 

4-4) ^(^.^). (4.10) 

Assuming that, as in the mean field case, the tip magnetization is proportional to the 
critical correlation function between u = and u = D ~ (^/C) 1 ^ (cf Figure 4.2), one 
obtains for a < 1 from (4.5) 



m(t) ~ exp 



(1-a) V C 



j.-v(l-a) 



1/a 



(4.11) 
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An essential singularity as encountered here has also been found in the six-vertex model 
[6] related to two-dimensional ice problems and to the roughening transition [59]. In this 
homogeneous system, however, all relevant quantities (free energy, correlation length and 
order parameter) vary in this way and the transition, which has been called of infinite 
order, is different from the present one. 



4.3. Ising model 

Calculations have been done to confirm in particular the behaviour (4.11) of the 
magnetization [56, 60]. For this the corner transfer matrix method was adapted to the 
parabolic geometry. Thus one works with the transfer matrix relating the spins along 
the upper and the lower boundary of the system. In the Hamiltonian limit this leads to 
an inhomogeneous quantum spin chain (Appendix B.2). If one deals with a lattice, the 
boundaries actually are step functions and the same then holds for the coefficients in the 
chain Hamiltonian (B3). Within a continuum interpolation they take the forms 

h n = (2n + fi) (X \ n = \(2n + v) a (4.12) 

so that h n and A n vary as n Q , reflecting the shape of the system. The parameters /it, v 
allow for differences in the total number of vertical and horizontal bonds at position n. For 
the case a = 1/2, /i = 0, u = 2, the chain Hamiltonian can be diagonalized with the help 
of special polynomials [56] . For other cases the problem was studied numerically and via a 
continuum approximation [60]. It was found that, for a < 1, the single particle eigenvalues 
uji vary as (21 — 1)/L 1_a for a finite critical system and as Z Q /£ 1_a near criticality, reflecting 
again the geometry. For a "paraboloid" the formula (B19) can then be used to obtain the 
magnetization at the tip. It reproduces the conformal prediction at criticality and gives 

m ~ £-(i-«)/2 exp L a £-(i-a)A*j (4.13) 

near the critical point. The exponential factor coincides with the expression (4.11) for 
v = 1. Formally, the exponential dependence appears, because the eigenvalues ui scale 
with a power of t in contrast to the case of a wedge where they vary logarithmically with 
t. 

Calculations for the case of free boundaries have not yet been done with this method. 
However, the problem was studied via Monte Carlo simulations and the same typical 
behaviour of m was found [57]. One expects a different prefactor in this case but the 
data did not allow its determination. Physically, the exponential vanishing of m can be 
understood as follows: a system with a < 1 is asymptotically narrower than any wedge 
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and so the magnetization close to the critical point should lie below an arbitrary power 
law. 

4.4. Other systems 

The previous considerations may be extended to systems displaying anisotropic critical 
behaviour. In this case the correlation lengths in two perpendicular directions diverge with 
different exponents v\\ and v±. Then rescalings have to be performed with different factors 
b± = b and b\\ = b z , where z = v\\/v±_ [61]. Choosing the symmetry axis of the figure along 
the direction with v\\ , the parameter C changes according to 



and marginal behaviour occurs for a = 1/z. Examples are provided by systems with 
uniaxial Lifschitz points [62], directed walks or directed polymers [63]. In the last two 
cases z = 2 and the borderline geometry is the normal parabola or, in three dimensions, 
the paraboloid. This situation was studied in [64, 65]. 



C = b 



zct — l 



c 



(4.14) 
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Figure 4.3. Exponent a for the directed self-avoiding walk as a function of the parabola 
parameter C. 



The simplest case is the directed self-avoiding walk in two dimensions. It is obtained by 
considering a usual one-dimensional random walk x = x(t) as a path in the (x, t)-plane. In 
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the continuum limit, the relative number P of walks between two points (x,t) and (x',t f ) 
follows from the diffusion equation via 



This is the analogue of the Laplace equation for the critical mean field correlation functions 
considered before. At x = ±Ct a a time-dependent boundary condition is posed. In 
particular, if paths reaching the surface are terminated there, a necessary condition to 
obtain equal weights for all the walks with N steps in the original problem, one comes 
back to the Dirichlet boundary condition [P = 0) considered before. 

In the marginal case a = 1/2, if the paths start at the tip, the solution for large t is [64] 



The confluent hypergeometric function F has to vanish at the boundary, from which the 
exponent a follows. It is shown in Figure 4.3. For C — > oo the usual Gaussian distribution 
with a = 1/2 is reobtained. As C is lowered, the exponent a increases monotonously, thus 
showing the expected non-universal behaviour. The decay of P with t thereby becomes 
faster in narrow systems, as found before for z = 1 in wedges and cones. Other exponents 
can be defined as usual. For example the survival probability, i.e. the relative number of 
walks reaching the time t (or having N steps in the discrete case) varies as t 7p_1 with a 
susceptibility exponent 7 P = 3/2 — a. 

For the case a < 1/2, the effect of the geometry is relevant and the form 



was deduced. It shows the same stretched exponential behaviour as the correlation 
functions in the previous sections. For a > 1/2, the boundary is in a region where P 
is exponentially small anyway and therefore does not affect the behaviour in agreement 
with the scaling contained in (4.14). The results for the directed polymer problem, obtained 
on a discrete lattice via transfer matrix methods, are very similar [65]. 

Finally, normal SAW's with parabolic boundaries have been studied [66]. It was found 
that, for a < 1, this originally isotropic system becomes anisotropic. A chain with N 
monomers then involves two radii of gyration: R\\ ~ N v w along the axis of the system and 
R± ~ N Ul - in the transverse direction. The chain configuration can be described using a 
blob picture [36] as shown in Figure 4.4. It results from the piling up of self-avoiding blobs 
inside the parabola. Within each blob the correlations are the same as in an unconfined 




(4.15) 




(4.16) 




(4.17) 
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Figure 4.5. Variation with a~ of the SAW exponents i>\\ and u± in two dimensions for a chain 
confined inside a generalized parabola. The points are the results of Monte-Carlo simulations 
on a square lattice. The surface geometry is relevant and induces an exponent anisotropy when 
a<l. 



chain. The radius R±_ is given by the transverse size of the system at a distance R\\ from 
the tip so that v± = av\\ . This implies a new anisotropic fixed point for which, according to 
(4.14), the geometry remains invariant under renormalization since z = 1/a. The evolution 
of the blob size along the axis of the system simply follows from the geometry and R\\ is 
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easily obtained as a function of N leading to 

"II = W 7 ( 4 - 18 ) 

11 a + -a) 

where z/ is the bulk exponent {y = 3/4 in d = 2). Monte-Carlo simulations on a square 
lattice (Figure 4.5) confirmed this picture and gave results compatible with (4.18). An 
accurate calculation of the correlation function like (4.16) has not yet been done. 



5. Extended defects at surfaces 

In the last sections, systems were treated where the geometry modifies the critical 
behaviour. We now turn to problems where a smoothly varying inhomogeneity in the 
couplings, decreasing from the surface, may influence the local critical properties. In real 
systems such an inhomogeneity could result from surface induced elastic deformations of 
the lattice. Thus one studies couplings which deviate from the bulk ones by an amount 

y being the distance from the surface. 

Under a scale transformation r' = r/b the coupling inhomogeneity, being an energy-like 
perturbation, transforms as AK'(y') = b 1 / u AK(y), therefore [67, 68] 

A' = b x l v -"A (5.2) 

where v is the critical exponent of the correlation length. Thus, for u > 1/v, A decreases 
and one expects the same local critical behaviour as at a free surface. On the other hand, 
in the marginal case (a; = 1/u) and in the relevant case (oj < 1/u) scaling theory predicts 
local properties different from those of a simple free surface. 



5.1. Ising model, general results 

An extended defect of the form given in (5.1) was first studied by Hilhorst and van 
Leeuwen [69] in the two-dimensional Ising model. They considered the problem on a 
triangular lattice at the critical point. The nearest neighbour couplings K\(y) parallel to 
the surface and the diagonal ones K 2 (y) were different from their bulk values as in (5.1) 
with amplitudes such that AK 1 (y)/AK 2 (y) = 2 sinh(2Ki)/ cosh(2if 2 )- 

The calculation of surface quantities was based on the repeated application of the star- 
triangle transformation [70] used also in calculations of the corner magnetization (Section 
3.2). From the evolution of the couplings one can deduce the surface magnetization m s 
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and the surface spin correlation function G\\ for any type of inhomogeneity in a numerical, 
iterative way. For smoothly varying perturbations as in (5.1), the asymptotic behaviour 
of the couplings after a large number of iterations and consequently the surface critical 
properties of the model can be determined exactly. 

The results for the critical behaviour [69, 71, 72] are in complete agreement with the 
scaling arguments in the previous section. In detail they are: 

i) For oj > l/u = 1, the perturbation is irrelevant, m s {t) ^t 1 / 2 and G\\(r) ~ r _1 as in 
the homogeneous, semi-infinite system. 

ii) In the marginal case, ui = 1, when A is smaller than a critical value A c > 0, the 
decay of Gu(r) is non-universal with 2x s = 1 — A/A C , while at A = A c where x s = the 
decay is logarithmic, G\\(r) ~ (logr) -1 . For strong enough enhancement of the couplings, 
A > A c , there is a spontaneous surface magnetization at the bulk critical temperature. 
As A approaches A c from above, it vanishes as (A — A c ) 1 / 2 . The spin correlations in the 
surface approach their limiting value m 2 s according to a power law with the non-universal 
exponent 2x s = A/A c — 1. 

iii) For uj< 1 the perturbation is relevant. For any A>0, there is a spontaneous surface 
magnetization at the bulk critical point which vanishes as A 1 /^ 1-1 ^ as A approaches zero. 
The correlation function G\\ has a stretched exponential form for any sign of A 

G|,(r) ~ exp [-(r/^i) 1 -] £|| ~ |A|-Vd— ). (5 . 3) 

The problem was also studied on the square lattice using Pfaffian methods [73, 74]. 
In this case the couplings parallel to the surface were chosen to be constant K\(y) = Ki, 
while the perpendicular ones were modified as in (5.1), writing 

AK 2 (y) = 4;, A = j sinh(2K 2 ). (5.4) 

For the marginal problem (u = 1) the complete set of surface exponents, both 
those defined at the critical point and those associated with the approach to criticality, 
was determined. Furthermore, including a surface magnetic field h 8 , the susceptibility 
exponents were also calculated. All exponents vary continuously with A. At the critical 
point the results are analogous to those obtained on the triangular lattice. The critical 
exponents have the same dependence on A/A c in both cases. This is a kind of "weak 
universality" as observed also in other models containing a marginal operator [75] . 

The critical behaviour is anomalous in the regime of spontaneous surface order A > A c . 
Then the scaling is anisotropic for T > T c and the local exponents are asymmetric, i.e. 
different for T < T c and T > T c , respectively. This unusual behaviour can be explained 
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[76] assuming that the local energy density is a "dangerous irrelevant variable" [77] for 
T > T c , whereas it is harmless in the region T <T C . 



5.2. Ising model boundary magnetization 
For the Ising model on the square lattice, the surface magnetization at zero surface 
field is easily obtained using the transfer matrix method (Appendix B.l) [78, 79, 33]. It has 
been calculated first for the homogeneous system and later for inhomogeneous ones. The 
quantity m s is deduced from the asymptotic value of the correlation function G\\. In the 
anisotropic limit from (Bll) it is found to be equal to the matrix-element of the operator 
erf between the ground state |0 > and the first excited state |1 > of the Hamiltonian in 
(B3) 

m s =< l|of |0 > (5.5) 

In the general case, this has to be multiplied by C\ = cosh K\ where K* is the dual coupling 
defined below (B2). This expression holds in the low-temperature phase where the state 
|1 > is degenerate with the ground state in the thermodynamic limit. 

Using the fermion techniques described in Appendix B.l one can express the matrix 
element by the surface component of the eigenvector $ s corresponding to the smallest 
excitation energy in (B8) 

< ljcrf jO >= $ s (l). (5.6) 

Since the vector <fr s must be normalized, its component $ s (l) is only non- vanishing in the 
thermodynamic limit if the eigenstate is localized near the surface. Such a state, with a 
localization length diverging at the critical point, is known to exist below T c [80, 78, 79]. 

Using (B8) and (B9) with h n = 1, the surface magnetization of a very anisotropic, 
semi-infinite system can be expressed in closed form 



m. 



CO I 

e n ^ 

1=1 n=l 



-1/2 

(5.7) 



Including a factor C\ and putting 

tanhifsH ,_ Q s 
A " " tanhK* • (5 ' 8) 

the expression also holds outside the anisotropic limit. Close to the critical point, $ s (n) 
extends a long way into the bulk, therefore its normalization and thus m s are determined 
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by its asymptotic behaviour for n ^> 1. For the couplings (5.4), in the extreme anisotropic 
limit, this gives 

An — A ^1 + ^n-^j (5.9) 

where A = 1 at the bulk critical point while A > 1 and A < 1 correspond to T < T c and 
T > T c , respectively. 

For the homogeneous model (A = 0) the sum in (5.7) is a geometric series. This gives 
m s = Ci(l — A~ 2 ) 1//2 which is the result first found by McCoy and Wu [3]. Near the critical 
point m s ~ t 1 ! 2 so that the magnetic surface exponent is f3 s = l/2. For to > 1 the products 
in (5.7) are convergent, thus for large n the behaviour of $ s (n) is the same as for A = 0. 
As a consequence (3 8 remains unchanged and the perturbation is irrelevant. When uj < 1, 
$ s (n) varies asymptotically as exp [— (A/2)n 1 ~ UJ ~\ corresponding to a localized state for 
A > 0. Then there is a spontaneous surface magnetization at the critical point which 
vanishes as A 1 /^ 1- ^ when A goes to zero, in accordance with results on the triangular 
lattice. For A < on the other hand, there is no surface order at the critical point and 
m s vanishes with an essential singularity according to 



m s ~ exp 



-a(A, uj)t' 



(5.10) 



In Figure 5.1 the temperature dependence of m s is shown for an isotropic bulk system and 
various values of A. 
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Figure 5.1. Boundary magnetization m s vs. temperature in the Hilhorst-van Leeuwen model. 
The full lines correspond to oi = l and various values of A, the dashed line to a>=0.25, A= — 1.5. 
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In the marginal case u> = 1, a closed form for m s can be obtained if (5.9) is changed 
into \ n = A(l - A/2U)- 1 . Then 



m, = 



FM - 1; A" 2 



-1/2 



(5.11) 



where F(a,b;c;z) denotes the hypergeometric function [81]. Expanding (5.11) in powers 
of £ = 1 — A -2 one obtains 



\/r(i - ^ 



t (l-A)/2 



A< 1 



(5.12) 



thus P s = (1 — A)/2. Expressions for A > 1 can be obtained in the same way. 



5.3. Conformal invariance for the marginal case 
In section (2.2) it was shown that the special conformal transformation (2.26) which 
leaves the surface invariant can be used to obtain information on the critical behaviour. 
For the system considered here the geometry is the same and the inhomogeneity (5.1) 
transforms into itself in the marginal case uj = 1/u [82]. This suggests that conformal 
invariance holds in the same way as for the homogeneous case. In particular one can then 
map the half-space onto a strip as in (A22) and investigate the specturm of the transfer 
matrix in this finite geometry. Under this mapping the inhomogeneity is transformed into 
a sinusoidal form [83] 

(5.13) 

where L is the width of the strip and -u, with < u < L, is the transverse coordinate. AK 
does not depend on the position along the strip. 

For the Ising model, a strip with couplings given by (5.13) was studied at the bulk 
critical point by two different methods. The case of a triangular lattice was treated in [83] 
by a numerical method based on the star-triangle transformation [84]. For the square 
lattice, the spectrum of the transfer matrix was obtained analytically in the extreme 
anisotropic limit [83]. Thus the Hamiltonian H as in Appendix B.l was investigated 
for couplings given by h n = 1, A n = 1 + AK(n)/2. The matrix equation (B8, B9) 
was transformed into a differential equation which is solved in terms of hypergeometric 
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Figure 5.2. Surface scaling dimensions for spin (<r) and energy density (e) in the Hilhorst-van 
Leeuwen model for uj=1 as a function of the strength parameter A. 

functions. In this way the lowest gaps, which are of order 1/L, were determined in a 
continuum approximation as in [85]. 

The spectrum found in the second approach displays the tower-like structure typical 
of conformally invariant systems. From the lowest gaps surface exponents are deduced, 
using (A23), which are in agreement with those found from the calculations for the half- 
space [73, 74, 86]. They are shown in Figure 5.2 as functions of the parameter A. For 
A > 1 the lowest gap vanishes like L~ A , i.e. faster than 1/L, in the large- L limit. Hence 
the smallest magnetic dimension vanishes. This is related to the fact that for A > 1 
there is magnetic surface order at the bulk critical point. These scaling dimensions are 
also in agreement with the results of the first-order perturbation expansion (Table C2 in 
Appendix C), where, according to (5.4) and Appendix B.3, the defect amplitude becomes 
g = A I '4 in the continuum limit. Similar conclusions are reached from the calculations 
on triangular lattices. These results strongly support the validity of conformal invariance 
for this particular inhomogeneous system. The correlation functions in the half-space 
therefore will have the same structure (2.34) as for a homogeneous system, but with a 
different scaling function. 

5.4. Related problems 
The Hilhorst-van Leeuwen model has been generalized recently by adding another 
marginal contribution with angular dependence to the inhomogeneity in (5.1). The 



36 



Inhomogeneous systems 



Hamiltonian spectrum in the associated strip was determined exactly, noticing the 
supersymmetric aspects of the eigenvalue problem [87]. 

For the most general marginal extended perturbation of the Ising model with the form 
Af(9)/p l / v in polar coordinates, the gap-exponent relation was shown to remain valid up 
to first order in the perturbation amplitude A [82]. 

In addition to the Ising model, the Gaussian model with a defect of the form (5.1) has 
also been studied [88] . This can be viewed as a mean field calculation which, however, is as 
complicated as the one for the Ising model. Using again the star-triangle transformation, 
the surface correlation function was determined as in Section 5.1. In the marginal case, 
which here corresponds to uj = 2, it decays with a non-universal exponent. 

For systems with anisotropic scaling as described in Section 4.4, the relation (5.2) has 
to be used with the exponent v in the direction of the inhomogeneity. For example, for 
a polymer directed along a surface and interacting with it, the monomer fugacity can be 
modified as in (5.1). Then one has u± = l/2 in the transverse direction and the marginal 
case corresponds to u = 2. For weak enough surface interaction (A < A c ), the polymer 
is free and can be characterized by critical exponents which depend on A. For stronger 
surface interaction, it becomes localized near the surface and the transition at A = A C is of 
infinite order [89]. 

The defects considered so far extended from the surface into the bulk. One may also 
consider perturbations which are confined to the surface but vary along it. For the Ising 
model the effect of an inhomogeneous field 

h(x,y) = ^-6(y) (5.14) 

was studied in [90]. Under rescaling, the amplitude transforms as 

A ' = b i-x s -u> A ( 515 ) 

where the scaling dimension of h s in (A6) was used. For the Ising model, x s = 1/2 (Table 
2.1), which leads to marginal behaviour when u> = 1/2. This also follows from Appendix 
C.2 with d e ff = x s . Mapping the system onto a strip then leads to a homogeneous 
boundary field of strength ^n/LA which shifts the levels of the Hamiltonian by an amount 
of order 1/L. The magnetic surface exponent determined via (A23) varies continuously 
between xi=x s = l/2 for A = and xi = 2 for A^oo. The latter case corresponds to fixed 
boundary conditions and is also realized when the perturbation is relevant. The correlation 
function decays towards a finite value in the strip and the exponent refers to the connected 
part. 



37 



Inhomogeneous systems 



A special form of inhomogeneity is present in quasiperiodic systems which interpolate 
between periodic and random systems and where the perturbation is spread over the whole 
volume. Results have recently been obtained for two-dimensional layered systems with 
quasiperiodic or aperiodic modulations. The relevance or irrelevance of the aperiodicity is 
connected to the strength of the fluctuations in the couplings [91] and the corresponding 
criterion is obtained [92] through a generalization of Harris's argumentation for random 
systems [93]. Exact results have been derived for the specific heat [91] and the surface 
magnetization [94, 92] for different aperiodic two-dimensional Ising models and their 
related quantum chains, as well as for directed polymers [92]. It was found that irrelevant 
aperiodicities indeed do not modify the critical behaviour, while for relevant perturbations 
the thermodynamic quantities show essential singularities at the critical point like in (5.10). 
As expected varying exponents are found in the marginal case. 



5.5. Scaling considerations for relevant inhomogeneities 
In the previous sections we have seen that the critical behaviour of different systems in 
the presence of relevant perturbations (u; < l/i/) is very similar: critical correlations decay 
in a stretched exponential form, while the thermodynamic quantities have an essential 
singularity at the critical point. In the following we show how a general scaling theory 
together with plausible assumptions can explain these observations. 

We start by considering the scaling behaviour of the surface magnetization 

m s (t,A) =b- x °m s {b 1/l, t,b l / v -"A) A>0 (5.16) 

where we included the relevant scaling field A, with the transformation law in (5.2). Putting 
b = £ in (5.16) one obtains m s in the form 

m s (t,A) = t (3 sf (^j l = \A\~ v ^- u ^ (5.17) 

where / is a characteristic length introduced by the inhomogeneity which stays finite at the 
critical point. If m s is to stay finite at bulk criticality the temperature dependence has to 
cancel. Then 

m s ~ a vx '^- v ^ A > 0. (5.18) 

For the Ising model this is the relation cited before in Section 5.1. 

The scaling behaviour of the parallel correlation function is given by 

G\\ (t, A, r) = b~ 2x 'G\\ (b 1/u t, b 1/u -"A, T -). (5.19) 
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Taking b = r, 

G|| (t, A, r) = r~ 2Xs g\\ (r^t, r 1 '^ A) (5.20) 

where the scaling function in the homogeneous case ( A = 0) behaves like g\\ (u, 0) ~ 
exp(— u v ). To obtain information on the case t = 0, A ^ consider the choice 10 = 0. Then 
one is dealing with an off-critical system where now A is playing the role of t. Therefore 
<7||(0, v) ~ exp(— v v ) in this case. Assuming that this also holds for ui ^ one arrives at 

G\\ (r) ~ r~ 2Xs exp [-(r/Z) 1 "^] . (5.21) 

for the critical correlations. This is the behaviour which was found in (5.3) for the Ising 
model [72] . A similar argument leads to the scaling form of the perpendicular correlation 
function at t = 

G ± (r) ~ r~ {x+x ^ exp [-(r/l) 1 '^] . (5.22) 

For reduced surface couplings [A < 0) the surface order, below the bulk critical 
temperature, is induced by the bulk magnetization at a distance D from the surface. The 
size of this surface region can be estimated by equating the thermal and inhomogeneity 
energy contributions, t ~ \A\/D UJ . One may then argue as in Section 4.2 and assume that 
the order parameter near the surface is proportional to the correlation function G\\ (r) with 
r = D ~ (\A\/tY^ so that 

m s (t,A) -exp [-alAI 1 /^- 1 ^] A < 0. (5.23) 

This is consistent with the Ising model result (5.10). 

The results (5.21-23) have the same functional form as (4.5) and (4.11) for the parabolic 
geometry with the correspondence a <-> lots, 1/C <-> \A\ U . This can be understood as 
follows: in both problems, a position-dependent, smoothly varying local length-scale can 
be defined near T c . For the parabola, it is proportional to the local width 

£(x) = Cx a (5.24) 

while for an extended defect, since AK(y) is a local temperature shift 

£(y) = &K(y)-». (5.25) 

For the calculation of a correlation function between boundary and bulk variables, one can 
imagine the system divided into sections, so that in each of them there is a decay with the 
local length or £(y). The complete correlation function is then obtained in the form 
of a product, or equivalently, as 

G , fol , w) ~ex P (-£!^). (5.26) 

Inserting for instance the form (5.25) then gives 

GiCvi.Ifc) ~ exp [-a (y\— - y\~^)] (5.27) 
which corresponds to (5.22). 
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6. Narrow line defects in the bulk 

In the last section the effects of extended defects at surfaces were treated. Narrow 
defects, where the couplings at the surface are modified only in a few layers do not change 
the critical behaviour in two dimensions. If the couplings are modified along a line, at 
a finite distance from the surface, a singularity in the surface correlation length appears 
below the bulk critical temperature [95]. This singular behaviour, however, occurs for a 
finite value of the correlation length and is not linked with any new fixed point. However, 
if defect lines are situated in the bulk, a rich local critical behaviour can be found. In this 
section we consider this situation with an emphasis on line and star defects. 

First some mean field and scaling considerations are presented, which show how the 
dimensionalities of system and defect, the bulk and the surface scaling dimensions may 
affect the local critical behaviour. This is followed by a discussion of line and star defects in 
the two-dimensional Ising model for which exact results have been obtained, either directly 
in the plane geometry or, using a conformal mapping, in the strip geometry. 



6. 1 . Mean field theory and scaling 
In mean field theory the problem can be discussed for arbitrary dimension d. Assume 
that the system contains a defect of dimension d* (a point, a line or a plane) so that the 
order parameter depends only on one coordinate r. Then the Ginzburg-Landau equation 
in the ordered phase reads in terms of reduced variables m = m/nib, f = r/£ [96, 97] 

i-d+d*+i d f ~d-d*-idm\ _ ^ 3 



r~ a+a +1 tt \f a - a ~ L =rr \=m 6 -m. (6.1) 
dr \ dr J 

The spontaneous magnetization takes the scaling form 

T 

m(r) = m b f(-) (6.2) 

as in (2.17) for the free surface. Equation (6.1) has to be supplemented by boundary 
conditions. If the defect interactions are weaker than bulk ones, the extrapolation length 
discussed below (2.3) is positive and one may look for the conditions under which a solution 
exists which satisfies Dirichlet boundary conditions at the defect. Since m < 1 near the 
defect, the local behaviour can be deduced from a linearized equation 

d 2 m d — d* — 1 dm „ .„ „ s 

-tt^ + : -7^ + m = (6.3) 

dr A r dr 
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which is of the Bessel type. An acceptable solution is obtained only when d — d* < 2. It 
takes the form f M J M (r) with 2/u, = 2 — d + d* and behaves like f 2m for small f so that, near 
the defect 

m(r) ~ r 2-d+d* t (3-d+d*)/2 d-d* <2 (6.4) 

The perturbation is then relevant: the defect changes the local critical behaviour and leads 
to a magnetization exponent /?/ = (3 — d + d*)/2. For d* = d — 1 one recovers the ordinary 
surface behaviour with $ = (3 S = 1. If d — d* > 2 the perturbation should be irrelevant 
and one expects bulk critical behaviour. This can be shown using scaling arguments as 
follows. The coupling perturbation which belongs to the defect subspace of dimension 
d* and couples to the bulk energy density e with dimension x = d — 1/u has a scaling 
dimension 

d*- x = d*-d+^ (6.5) 

As a result the perturbation is relevant (irrelevant) when d — d* < 1/u (> 1/v). With the 
mean field value v = 1/2, one finds the condition cited above. 

If d* = d— 1, the defect divides the system into two parts. Then if v < 1 and the 
local interactions are weaker than in the bulk, according to (6.5) the defect is expected 
to behave in the same way as a cut. Thus ordinary surface critical behaviour will result 
[15] provided the corresponding surface fixed point remains stable against a weak coupling 
between the surfaces of the two half-spaces. The condition for this can be found in the 
following way: starting from two free surfaces the perturbation involves the product of 
two surface magnetization operators with scaling dimension x s and the dimension of the 
coupling is given by 

d-l-2x s = ^t (6.6) 

where 7 S is the local susceptibility exponent at the ordinary surface transition (A9). The 
system with the defect behaves in the same way as one with a free surface if this exponent 
is negative and v < 1 [98]. Non- universal local critical behaviour is expected when the 
perturbation is marginal, i.e. when v = 1 and 7 S = 0. These scaling results have been 
confirmed in a series of calculations on the O(N) model in the limit N — > oo and using 
either e- or 1/A-expansions [99-107]. 

6.2. Ising model with a defect line 
According to the scaling arguments of the preceding section a defect line in the two- 
dimensional Ising model, which has exponents v = \ and 7 S = at the ordinary transition, 
is expected to lead to continuously varying local exponents. This problem was first 
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Figure 6.1. Ising square lattice with two kinds of defect lines (a) chain defect, (b) ladder defect. 

investigated by Bariev [108] who deduced the local magnetization, as a function of the 
perturbation strength and distance to the defect, from the asymptotic behaviour of the two- 
spin correlation function below the critical temperature. This was followed by a detailed 
study of the two-spin correlation function by McCoy and Perk [109]. 

In the square lattice Ising model one can distinguish two types of perturbation as shown 
in Figure 6.1. The chain perturbation has modified couplings K[ parallel to the defect line 
whereas, for the ladder defect, perturbed couplings K' 2 are in the perpendicular direction. 
The local magnetization was found to vary as [108] 

< a(y) >~t pl y pl -P (6.7) 

when the distance y to the defect line is much smaller than the bulk correlation length £. 
In this expression (3 = 1/8 is the bulk magnetization exponent and (3i the continuously 
varying local magnetization exponent. The spatial variation in (6.7) simply follows because 
< a(y) > has the scaling form (6.2) with v = 1. 

The local magnetization exponent varies with the type and strength of the defect 
according to [108,109] 

A = 4, arctan 2 «i = e" 2 ^"^ = chain defect ( 6 - 8 ) 

A = \ arctan 2 ^^" 1 ) k 2 = ^hif 2 ladder defect (6.9) 
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(a) (b) 



Figure 6.2. Local magnetization exponent /3; vs. defect strength in the Bariev model for (a) 
chain defect and (b) ladder defect. 

where the bulk couplings have their critical values and the asterix denotes dual variables 
(see below (B2)). These exponents are shown in Figure 6.2. 

The chain defect exponent decreases continuously from (3\ = 1/2 when K[ — > — oo to 
Pi = when K[ — > +oo (Figure 6.2a). In the first limit the spins along the defect are 
frozen into antiparallel configurations and nearby spins couple to a vanishing total spin. 
The local exponent then takes the free surface value. In the other limit the vanishing 
asymptotic value is linked to the onset of local order at the bulk critical point. When 
K[ = one can sum out the inner spin and the chain defect effectively becomes a ladder 
with strength given by tanhi^ = tanh 2 X 2 . Then (6.8) and (6.9) give identical results 
since sinh 2Ki sinh 2K 2 = 1 on the critical line. 

For a ladder defect the local exponent is invariant under the change K' 2 — > —K' 2 since 
the original defect coupling can be restored through an appropriate spin reversal in one 
half of the system. In the positive sector (Figure 6.2b) f3i decreases from the free surface 
value when K 2 = to a nonvanishing limiting value when K' 2 — > +00. Then the ladder 
is also a chain defect with K[ = 2K\ and, with the appropriate values of the perturbed 
couplings, the two formulae give identical results. 

In their study of the spin-spin correlation function parallel to the chain defect McCoy 
and Perk [109] showed that the correlation length exponent keeps its bulk value v = 1 
although the amplitude of £y may depend on the local interaction strength. At the critical 
point the parallel correlations decay with an exponent rf\\ = 2xi where xi = (3i as expected 
from scaling, compare (A7). 
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The incremental specific heat introduced by a chain defect in the square lattice Ising 
model was obtained by Fisher and Ferdinand [110]. It has a universal exponent a\ = 1 
corresponding to the expression (A2) with d replaced by d*. 

The energy density correlations were calculated exactly in the scaling region in [111]. 
Although the correlation function displays much structure, the decay exponents are 
universal, keeping their bulk values. The same conclusion was reached previously in 
[112, 113] using the techniques of operator algebra. Although the line defect introduces 
a marginal operator which is the energy density, this operator has to keep its universal 
scaling dimension x = 1 in order to obtain a continuous variation of (3i with the defect 
strength: Suppose the dimension of the energy density operator to be changed to x(AK) 
through the introduction of a defect AK. Then a further change of the defect strength 
would no longer constitute a marginal perturbation since its scaling dimension, given as 
in (6.5) by 1— x(AK), would be nonvanishing. 

Burkhardt and Choi [114] obtained the form of the critical n-point correlation function 
for the energy density by summing a perturbation series in the defect strength. The 
complete finite-size behaviour of the transfer matrix spectrum below T c was investigated 
in [115]. 

The Hamiltonian limit of the line defect problem has been studied through real-space 
renormalization [116] and low-temperature expansion [117]. In this limit one can see a 
close connection between the line defect problem and the X-ray problem [118] in which 
correlation functions with continuously varying exponents also occur. Using this analogy 
[119] one obtains the decay exponents in agreement with (6.8) and (6.9) where, in the 
Hamiltonian limit (Appendix B.l), 

K'* K' 

The size dependence of the lowest gap in the ordered phase has been investigated 
through numerical and analytical calculations [120-123]. 

A one-dimensional quantum hard-dimer model, belonging to the Ising universality 
class, has also been studied [124]. The novel feature in this case is the jump, with increasing 
defect strength, from non-universal to ordinary surface behaviour for a critical value of the 
defect coupling. 

Nigthingale and Blote [125] used finite-size scaling on a strip to study line defects in 
2d models belonging to the Ising universality class as well as in the (/-state Potts model. 
Varying magnetic exponents are obtained in the Ising universality whereas the Potts model 
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displays bulk behaviour for q = 1/2 (irrelevant perturbation) and strong crossover effects 
for q — 3 (relevant perturbation). 

The Potts model with a defect line has also been studied through exact renormalization 
on hierarchical lattices [126]. In the relevant case, with strengthened couplings, the defect 
remains ordered at the bulk critical point and the system displays a local first order 
transition [76]. 

6.3. Conformal invar iance, star-like defects 
For a system with a line defect the symmetries usually considered to be necessary 
for conformal invariance (Appendix A. 2) are partially broken. This is a similar situation 
as for a surface. Nevertheless the gap-exponent relations and other spectral properties 
typical of conformally invariant systems still survive in the marginal Ising case. This 
was first conjectured on the basis of a numerical study of the classical system mapped 
onto a strip [127]. The Hamiltonian limit was then studied for larger sizes using fermion 
techniques [128]. The form of the whole spectrum was conjectured in [129] and exactly 
calculated in [130]. 



3 2 



A. 




z- plane w-plane 

Figure 6.3. A star of defect lines in the z-plane and the corresponding configuration in the 
periodic strip of the w-plane. 

Under the logarithmic mapping (A17) a single infinite defect line in the plane is 
transformed into a pair of parallel and equidistant lines along a strip with periodic 
boundary conditions [127]. The same mapping can also be used for star-like defects in 
the plane [131]. A star with rid semi-infinite linear branches at the polar angles 6j = 2iv5j 
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(j = 1, . . .rid) is then mapped onto rid parallel defect lines in a periodic strip of width L 
(Figure 6.3). With ladder defects, in the extreme anisotropic limit, the Hamiltonian (B3) 
becomes 

1 . L ~ 1 _ nd . 

« = - 2 E \-< + + 2&- ( 6 - n ) 

n=l j=l 

where the defect strength kj is the ratio of perturbed to unperturbed couplings like in 
the second relation (6.10). The defect positions Lj scale with the size of the system as 
= Lj — Lj-i. 

This Hamiltonian has been diagonalized in [130] using fermion techniques (Appendix 
B.l). Between the defects the eigenvectors of the excitation matrix have a plane-wave 
form. The values of the wave-vectors follow from the 2nd boundary conditions at the 
defects. The conformal properties are preserved only when the Sj are commensurate. If to 
is the smallest common denominator of the S^s, either m = rid for equidistant defects or 
m is the number of defects in a corresponding equidistant configuration which is obtained 
by adding m — rid lines with Kj = 1 (Figure 6.3). Then one gets 7(777,) fermion species 
in the diagonal Hamiltonian. In the plane this is the number of lines crossing (m even, 
7 (to) = to/2) or meeting (m odd, 7(m) = to) on the star defect in the corresponding 
equidistant configuration. 

The 0(L _1 ) part of the Hamiltonian can be written in diagonal form as 



i=\ \r=0 



{r+\ -Ai)a+a ir + {r + \ +A l )(3+(3 l 



1 

12 



I -6A? 
2 



(6.12) 

where af r (oti r ) and (3f r (Ar) are fermion creation (anihilation) operators. The shifts Aj 
depend on the k/s and also involve the parity eigenvalues p = ±1 due to the periodic 
boundary conditions (see (B6)). When L is even the magnetization and energy sectors 
of the original Hamiltonian are in correspondence with the odd states of Ti-i and the 
even states of H+i, repectively. Due to the shifts Aj the magnetization sector contains an 
infinite number of conformal towers. 

In the case of a single defect line in the plane the two equidistant lines on the strip 
have the same strength k. Then m = 2, 7(771) = 1 and one gets [130] 

A(k) = 1 - f arctan^- 1 ) p = -1 (6.13) 

A(«) =0 p = +1 (6.14) 

The lowest gap gives the dimension of the local magnetization via (A20) and as in (A7) 

Pi = 1[A(k) - l] 2 = A arctan 2 ^- 1 ) (6.15) 

2 TV 2 
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in agreement with the direct calculation in the plane (6.9). When there are several defects, 
Pi refers to a star and is something like a generalized corner exponent depending on 
all the defect strengths. In the energy sector the spectrum remains unaffected by the 
defects and the scaling dimensions are universal. According to (6.12) when m > 2 the gap 
normalization differs by a factor 7(771) from the usual one for periodic boundary conditions. 
This point is discussed in [132]. The spectrum-generating algebra has been obtained in [133] 
for the single line problem and in [131, 134] for general values of m. Details about this 
aspect can be found in the review [135]. 

In addition to the defect problems described so far, the following situations also have 
been studied. For a 3-state quantum chain with a defect, corresponding to a Potts model, 
the structure of the spectrum was studied numerically in [136]. A defect where the Pauli 
matrix a x : , in one coupling term of the periodic quantum Ising chain, was replaced by 
a general Hermitian 2x2 matrix was studied in [137]. Exact results could be obtained 
when the global Z 2 symmetry of the quantum chain is preserved. Then the critical and 
conformal properties are the same as with an ordinary defect. The same occurs for the 
Ising quantum chain with staggered 3-spin interactions and the Ashkin- Teller model with 
generalized defects [138, 139]. Finally, the Ising quantum chain associated with star-like 
defects was generalized in [140] introducing modified couplings over an extended region. 
In order to keep the system critical, the same defect strength was taken for the a z and 
a x parts in the strip Hamiltonian (B3), which corresponds to a constant change in the 
anisotropy inside this region. Conformal invariance is preserved, as above for stars, for 
commensurate configurations. For an arbitrary varying defect amplitude, the equidistant 
level structure of the spectrum is maintained at high energies only. 



7. Extended line defects in the bulk 

For defects in the bulk, interesting situations do not only arise if one deals with narrow 
ones as in Section 6 but also if they are extended as (5.1). This type of problem was first 
introduced and investigated by Bariev [141-144]. For such extended defects originating 
from a core of dimension d* in a rf-dimensional system, scaling considerations give two 
separate relevance- irrelevance criteria. One is associated with the extended part of the 
defect and is the same as for the surface case. Thus the perturbation is relevant, irrelevant 
or marginal for a; < 1/v, u> > 1/v and to = 1/z/, respectively. The other criterion, related 
to the core of the defect, is the condition on the sign of d* — d + 1/v found below (6.5). 
The most interesting situation occurs when both conditions predict marginal behaviour. 
In this case, although the exponents are varying with the defect strength A, there are 
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discontinuities at A = [145]. Connected with that, a perturbation expansion for the 
exponents is divergent [141]. This feature is outlined in Appendix C. 



7.1. Ising model, conformal results 



In the two-dimensional Ising model an extended line defect with decay exponent uj = 1 
satisfies both marginality conditions since v = 1. An exact treatment of this problem, 
however, is more difficult than for the extended surface defect in Section 5, since most of 
the analytical methods used there cannot be adapted easily to internal defects. The only 
exception is the method based on conformal mapping. 

The problem was treated [145] for a square lattice as in Figure 6.1 with a defect centered 
at y = and couplings K2(y) varying as 



where A is defined in (5.4) and the bulk couplings Ki, K 2 have the critical values. When 
K 2 (0) = 0, the system separates into two uncoupled semi-infinite ones, with a surface 
inhomogeneity as studied in Section 5. 



K 2 (y) = K 2 + 



A 



(7.1) 



y 



X 




- 1 







A 



Figure 7.1. Local scaling dimensions for spin (<r) and energy density (e) for an extended defect 
with u=l in the bulk of an Ising model as a function of the strength A. 
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Using the logarithmic transformation (A17) the problem is mapped onto a strip. The 
inhomogeneity thereby becomes periodic with period L/2 and is given by 

(7.2) 

In each of the intervals, up to a factor of two, this has the same form as (5.13). 

As in Section 5.3 the lowest gaps were determined in a continuum approximation. Due 
to the different boundary conditions, the eigenvalues and hence, the local critical exponents 
are different from those in the surface defect problem. The scaling dimensions are shown 
in Figure 7.1 as functions of A. 

For enhanced local couplings (A>0) the energy of the lowest excited state, Ei—E ~ 
L~ 1 ~ 2A , vanishes faster than 1/L so that the scaling dimension of the magnetization is 
x\ =0 and the defect region remains ordered at the bulk critical point. On the other hand 
for reduced couplings (A<0) there is no order at T = T C and the local exponents are linear 
functions of A. The magnetic exponent is the same as for two disconnected semi-infinite 
systems (see Figure 5.2). At A = it is discontinuous. As a consequence a finite-order 
perturbation expansion for the gaps around the homogeneity point A = is not possible. 
The expansion for the first gap starts as [145] 

9-77-/1 A \ 

E 1 -E = T ^---logL + ...j (7.3) 

i.e. the first-order correction term is divergent, but the singular contributions sum up to 
a regular term in infinite order. This observation is in agreement with the results of the 
perturbation theory for the system in the original geometry (Appendix C). 
The critical exponent for the energy operator given by 

xi = l + \A\ (7.4) 

is continuous at A = where it reaches its bulk value. It can be determined independently, 
working in the plane and using finite-size scaling for the matrix-element of the local energy 
operator as in (A21). The coincidence of the two results supports the assumption that the 
system is conformally invariant at the critical point. 



AK(u) = A 



L (2itu 
— sin — — 
2tt \ L 
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8. Radially extended defects 



The systems studied so far had essentially all a layered structure. As a final example, we 
now examine the case where the inhomogeneity extends radially. The centre may be located 
either in the bulk or on a free surface. Only marginal temperature-like perturbations 
corresponding to a change in the interaction strength are considered. After some scaling 
considerations, exact results for the local magnetization in the two-dimensional Ising model 
are presented. This is supplemented by a discussion of the conformal aspects. 



According to the scaling arguments of Section 6.1, a short-range perturbation 
introduced by a point defect with d* = into a (/-dimensional system has a scaling dimension 
—d + 1/u which follows from (6.5). Since 1/u < d at a second order phase transition this 
type of perturbation is always irrelevant. This can be verified explicitly for the Ising 
model through a calculation of the local magnetization which may be expressed in terms 
of averages for the unperturbed system [146]. Following the method of Section 2.2 and 
using the rotational symmetry one can also obtain the critical correlation functions. They 
are completely determined and identical to the unperturbed ones. 

In the case of an extended defect and a perturbation amplitude decaying as a power 
of the distance from the centre, the relevance-irrelevance criterion is the same as for 
perturbations decaying from an interior line defect or a surface (Sections 5, 7 and 
Appendix C). For a temperature-like perturbation, according to Equations (5.2) and (C7), 
the amplitude has a scaling dimension 1/u — u where u> is the decay exponent of the 
perturbation. In the following we consider the marginal case where u> = 1/u. 

8.2. Ising model 

For the Ising model containing an extended defect with w = 1 exact results were 
obtained using the corner transfer matrix method [147]. 

The section of a square lattice forming the corner transfer matrix is shown in Figure 
Bib. For the present problem one chooses anisotropic couplings Ki and 



where n is the row index increasing from the centre. The Hamiltonian limit leads to a study 
of an inhomogeneous quantum chain (B3) with couplings h n = 2n, X n = X(2n + a + 1). 
As it stands, the corner transfer matrix describes a segment of an anisotropic system with 



8. 1 . Scaling considerations 




(8.1) 
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Figure 8.1. Local magnetization exponent /3; vs. defect strength for a radially extended defect 
(6) in the bulk and (s) on the surface of an isotropic Ising model. Here a is equal to -wot. 

opening angle 90°. However, as explained in Appendix B.3, it also describes an isotropic 
system for which the opening angle is infinitesimal. Therefore an isotropic system with 
rotational symmetry and opening angle 9 can be described by the transfer matrix 

T{9) = exp (- 1 enj (8.2) 

According to (B19) the magnetization at the centre involves the single particle 
excitations uj\ of the quantum chain. The critical excitation spectrum can be obtained 
by taking the continuum limit. Then the very centre of the system has to be excluded 
[148] . Taking a cut-off at r and an outer radius R, the critical excitations are given by 



where the z\ are solutions of a, \n(R/r) tanz = — 2z. An additional bound-state-like solution 
occurs when a < 0. The local magnetization exponent f3i then follows from the finite-size 
behaviour of the centre magnetization which is analogous to (A8). This gives 

n a 2 f 00 , sinh 2 w . t . 

2 J sinh(7r | a \ cosh-u) 

A bound-state contribution | a | /2 must be added when a < 0. The variation of the 
magnetization exponent with the perturbation strength is shown in Figure 8.1. The 
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exponent decreases smoothly when the defect amplitude is increased but, contrary to the 
case of surface or extended line defects, it vanishes only asymptotically. Thus there is no 
local order at the bulk critical point for any finite value of the defect strength. 



8.3. Conformal considerations 
The previous Ising results can also be deduced from a mapping of the perturbed critical 
system onto a strip [149, 150]. This mapping can done for an arbitrary conformally 
invariant two-dimensional system containing an extended perturbation with rotational 
symmetry. Then, in the continuum limit, 

-I3U = -(3Hc + g J r _h V(r, 6)rdrd6 (8.5) 

where ip(r,6) is some local field with scaling dimension In the following, as before, 
we assume that the decay exponent takes its marginal value ui = 2 — x^. Under the 
transformation (A17) the dilatation factor just compensates the decay of the perturbation 
in the marginal case. Therefore, in the periodic strip of width L, one obtains an 
homogeneous L-dependent deviation from criticality 

A = 9 (8.6) 

for the couplings which are conjugate to ifj. One should note that such a homogeneous 
deviation from the bulk critical couplings leads to a massive excitation spectrum as in 
(8.3) and, contrary to the case of extended line defects treated in Sections 5 and 7, the 
tower-like, equidistant-level structure of the spectrum is lost here. 

Let 0, with bulk scaling dimension x, be either the energy or the magnetization density. 
Its local scaling dimension at the defect, xi, can be deduced from the corresponding gaps for 
the strip through (A20). According to finite-size scaling, the perturbed gap G<f, transforms 
as 

G^A, L) = L-'G^L^A, 1) = L- X U [(2n) 2 -^g] (8.7) 
where is a universal scaling function [151]. This leads to the local exponent 

xi = [W 2 - X *9] (8-8) 

It depends on the perturbation amplitude g, as expected for a marginal perturbation, 
but does it in an universal way [150]. Different systems belonging to the same class of 
universality will show the same dependence on the perturbation amplitude. 
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For the Ising model, the perturbed gaps (or gap scaling functions) can be calculated 
either for an isotropic strip or using the Hamiltonian limit [152-155]. In the latter case the 
defect amplitude has to be rescaled as discussed in Appendix B.3. On a square lattice with 
modified interactions K(r) = K c + A/r one recovers the local magnetic exponent given in 
(8.4) with a = 8A [149]. For small values of A, it has the expansion 

(3i = xi = I - 2A + 161n2 A 2 + 0(A S ) (8.9) 
8 

To first order this is identical to the perturbation result in Table C2 with g = 2A. The 
factor of two in the continuum parameter g reflects the two bonds per unit cell on the square 
lattice. Due to the dual symmetry of the Ising model [156] the exponent of the local energy 
density has an expansion involving only even powers of the perturbation amplitude and 
reads [150] 

xi = 1 + 32A 2 + 0(A 4 ) (8.10) 

The same procedure applies with minor changes when the defect centre is located on a free 
surface. This was also studied for the Ising model [149]. To first order the result is again 
in agreement with the perturbation approach (Table C2). The corresponding magnetic 
exponent is also shown in Figure 8.1. 

9. Conclusion 

In this review we have presented a number of systems and situations where the 
universality which is normally found at continuous phase transitions can be absent. 
Superficially, the examples look quite different, but in each case the tendency towards 
order is modified locally in a particular systematic way. This is the basis for the common 
features one observes. 

It is not difficult to understand the results qualitatively. For example, if order is 
favoured this leads to a steeper magnetization curve and will tend to reduce the magnetic 
exponent. This explains the general form and the similarity of the graphs in the various 
marginal cases. There are, of course, some differences in details. The close parallel 
between effects from the geometry and from the modification of interactions inside a 
system is nevertheless remarkable. A similar relation is found if one compares homogeneous 
integrable systems at and off the critical point [35, 157, 158, 48]. In this case the corner 
transfer matrix can be used as a link. For the present systems the direct connection is less 
obvious. 

In two dimensions conformal mappings play an important role in the studies. The 
results show that they can be used even beyond their obvious domain. The gap-exponent 
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relation also holds in cases like in Section 8 where no tower of equidistant levels exists. 
This is an aspect which might deserve further study. The same is true for the situation 
in three dimensions where only a few, essentially mean field results have been obtained so 
far. 

The last point is also important since it should be possible to measure some of the 
effects. The most obvious experiment would be to look for the local order in systems 
with wedge-like, conical or parabolic shape. This is not easy since the local order will be 
smaller than inside the bulk, and already for planar surfaces such measurements are rare. 
Nevertheless they should be feasible with proper local probes and would add an interesting 
aspect to the general picture of critical phenomena. 
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Appendix A. Scaling and conformal invariance 

Near a second order phase transition the singular parts of thermodynamic quantities 
vary as powers of the parameters (scaling fields) t ~| T — T c |, h, etc. measuring the 
deviation from the critical point. This type of behaviour is linked to the self-similarity 
of critical fluctuations inside the correlation volume £ d where £ ~ t~ v is the correlation 
length. The system is then covariant under a global change of the length scale and singular 
quantities are homogeneous functions of their arguments. These properties form the basis 
of the scaling hypothesis [159] which lead to scaling laws relating critical exponents to 
a small number of fundamental ones. The scaling behaviour was also the source of the 
renormalization group ideas [1], allowing a calculation of the fundamental exponents. 

At the critical point, £ diverges and the system becomes scale invariant. This, together 
with more usual symmetries (rotation, translation), leads to the invariance under conformal 
transformation, first used in field theory [160] and later fully exploited in two-dimensional 
critical systems where it allows an exact determination of the critical exponents and much 
more [161, 16, 162] . 



A.l. Scaling and critical exponents 
According to scaling theory, when the lengths are rescaled by a factor b > 1, i.e. 
when r — > r/6, the scaling fields (like t and h) are changed by a factor b d ~ x where d is 
the dimension of the system and x the scaling dimension of conjuguate quantities (e.g. 
the energy density for t, the magnetization density for h). When d — x > (< 0) the 
corresponding scaling field grows (decreases) under rescaling. Such a field is said to be 
relevant (irrelevant) whereas it is marginal when x = d. The system becomes invariant 
under rescaling only when the relevant scaling fields vanish which corresponds to the critical 
point. 

Since irrelevant variables finally vanish under rescaling, only relevant and marginal 
scaling fields influence the critical properties, the marginal ones generally leading to varying 
exponents. 

Near the critical point, the singular part of the free energy density is a homogeneous 
function of the scaling fields and transforms according to 

f h (t, h, = b~ d f b b d ~*h, £j . (Al) 

In (Al) we kept only relevant variables and included the inverse of the linear size of the 
system as a new relevant scaling field, since at its bulk critical point a system is truly 
critical only in the thermodynamic limit, i.e. when 1/L = 0. 
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The critical behaviour of conjuguate quantities and their derivatives can be deduced 
from (Al) taking advantage of the arbitrariness of the dilatation factor b. The 
corresponding exponents are all related to the basic ones x, v through scaling laws and 
involve d. 

For example, the specific heat exponent a can be obtained from the second t-derivative 
of both sides of (Al) at /i = 0, 1/L = and taking b = t~ v with the result 

a = 2 - dv (Al) 

The magnetization m = —dfb/dh transforms as 

m (t, h, = b~ x m (b 1/v t, b d ~ x h, (A3) 

and the exponent /?, defined through m ~ t@ : can be obtained from (A3) taking h = 0, 
1/L = and b = t~ v as before, so that 

= vx. (A4) 

At the bulk critical point t = 0, h = 0, critical singularities are suppressed by a finite L 
value. Finite-size scaling exploits the way they develop when L — * oo in order to determine 
the critical exponents. For example, from (Al) with b = L, the free energy density behaves 
as L~ d at criticality whereas the magnetization density in (A3) gives 

m(L) ~ L~ x . (Ah) 

When the system is limited in space through a surface, besides the bulk term f b in 
(Al), new local contributions (surface, corner, etc.) to the free energy density appear, with 
singular parts behaving as above for the bulk. For example, the contribution of a surface 
with dimension d — 1 transforms according to [4] 



/. (t,h s , i) = b~^f s ^/%b d -^h s , ^ . 



(A6) 



where a surface magnetic field h s has been included. Repeating the previous argument, 
the scaling relations (A4-5) remain valid with the appropriate surface exponents replacing 
bulk ones, 

P* = vx s , (A7) 
m s (L)~L- Xs . (A8) 
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The surface susceptibility exponent 7 S defined through d 2 f s /dh 2 s ~ t~ la is obtained as 

la = v{d-l-2x a ). (A9) 

The bulk two-point correlation function for an operator ifj is obtained taking a functional 
derivative of the free-energy with respect to the appropriate position-dependent field 

G(r, ( ) S <*(0)*(r)>= 5i g 5J . (A10) 

If e.g. ip(r) stands for the magnetization operator, then its average < ip(r) >= m(r) gives 
the local magnetization. The transformation law of the two-point function then follows 
from (A10) and (Al) 

G(r, t) = b~ 2x G b 1/u tj (All) 

where x is the bulk scaling dimension of ip. At the critical point, t = 0, taking b = r, the 
power- law decay of correlations follows 

G(r,f = 0) = A {A12) 

Surface correlations can be investigated in the same way. The correlation function with 
two points at the surface behaves as 

G\\(r,t = 0) ~r- 2Xs (A1S) 

while the perpendicular correlations decay like r -( x + x s) ( see Section 2.2). 



A. 2. Conformal invariance 

Covariance under conformal transformations is expected to hold at the critical point of 
systems with short range interactions, which possess translational and rotational symmetry 
and are invariant under uniform scaling. When some of the above symmetries are broken by 
a marginal perturbation (e.g. for inhomogeneous systems) some of the properties associated 
with conformally invariant systems, like the gap-exponent relation, can be preserved as 
observed in specific examples. 

A conformal transformation r — > r'(r) can be seen as a generalization of uniform 
scaling, where the structure of the lattice is locally preserved, but the rescaling factor 
6(r) becomes a smooth function of the position. It follows from the Jacobian of the 
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transformation as b(r)~ d = det(dr'/dr). Since local fields transform as h(r) — > h'{r') = 
b{r) d ~ x h{r) the two-point function in (A10) transforms like 

< ^(ri)ip(r 2 ) >= b( ri )- x b(r 2 )- x < V>K)?/>K) > (AU) 

under a conformal transformation. This is a straightforward generalization of (All) with 
t = 0. Similarly the transformation law for an operator profile < ip{r) > is obtained as 

< V(r) >= b(r)~ x < i;(r') > . (A15) 

The conformal group for dimensions larger than two is finite-dimensional and contains 
rotations, uniform dilatations, translations and inversions. The special conformal 
transformation with an arbitrary translation a, 

^ = ^ + a, (A16) 

is constructed from the two last ones. Invariance under this transformation fixes the form 
of three-point functions [160] like scaling does with the two-point functions in (A12) and 
(A13). If a in (A16) is an infinitesimal surface vector, then one can use this transformation 
to find restrictions on the form of surface correlations, as shown in Section 2.2. 

The method of conformal invariance is especially powerful in two dimensions where the 
conformal group, being isomorphic with the group of complex analytic functions, becomes 
infinite-dimensional and strongly restricts the possible values of critical exponents for a 
broad class of systems [161, 163, 16, 162]. 

In two dimensions one may also use a complex mapping w(z) to go from one geometry 
to another [164]. When some critical correlation is known in the first geometry, it can be 
transformed into the second one. The local dilatation factor is then b(z) = \dwfdz\~ 1 . 

Two basic geometries are connected by the logarithmic transformation 

w = ^- lnz (A17) 

which maps the infinite z-plane onto a periodic strip of width L in the w-plane. 
Transforming the correlation function in (A12) according to (A14), one obtains 

(2tv/L) 2x 

< uiM«2, v 2 ) >= t- -2^7 ^7 ( Al8 ) 

[2cosh^(-ui -u 2 ) -2cos^(t>i - v 2 )\ 
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in the strip geometry where u measures the distance along the strip and < v < L denotes 
the transverse periodic coordinate. An expansion of the r.h.s. of (A18) for u\ > u 2 leads 
to 



2TT °° 

< ip(ui,vi)ip(u2, v 2 ) >=( — ) 2x ^2 amCWexp 



2tt 

■ — {x + m + m)(ui - 112) 

LJ 



x exp 



-y-(m - m)(vi - v 2 ) 

1j 



(A19) 



where a m = F(x + m)/T(x)ml. The strip correlation function can be also determined 
through a direct calculation using the transfer matrix method as described in Appendix B. 
Then, comparing (A19) with expressions like (B13), one reaches the following conclusions: 

i) The eigenstates of the Hamiltonian in (Bl) are labelled by pairs of integers (m, m), so 
that the spectrum exhibits a tower-like structure with energies E + (x + m + m)2ii /L 
and momenta given by (m — m) 2tt/L. 

ii) The gap between the ground state and the lowest excited state for which the matrix 
element of ip is non-vanishing, is related to the scaling dimension x via 

2n 

E 1 -E = —x. (A20) 
Lj 

This is the gap-exponent relation. 

iii) Finally, the matrix element itself is given as 

2tt x 3 



<1H0>=^ T J (A21) 

For the magnetization, this is in agreement with the finite-size scaling result (A5). 

To study surface correlations, the half-plane (y > 0) can be mapped onto a strip with 
free boundaries through the conformal transformation 

w=-\nz. (A22) 

7T 

Now the eigenstates of the Hamiltonian in the strip are labelled by an integer m = 0, 1, . . ., 
the spectrum has a tower-like structure with energy eigenvalues E Q + (x s + 111)71 / L and the 
gap-exponent relation reads 

E 1 -E Q = ^x s . (A23) 
Finally, the matrix element for ip near the boundary is given by 

<1|V> S |0>= (£j Xa . (A24) 
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Appendix B. Transfer Matrices 

In the transfer matrix method [156, 165, 166] a system of spins or other classical 
variables with short-range interactions is built from identical smaller units. This technique 
can also be used for the inhomogeneous systems considered here. In two dimensions one can 
choose the basic units as shown in Figure Bl. The corresponding transfer matrices T are 
the partition functions of either two consecutive rows (case a) or a whole angular segment 
(case b). Their matrix character results from their dependence on the two sets {<r}, {a'} of 
boundary spins (full points). The partition function Z of the whole system is then obtained 
by taking the product of an appropriate number of T's. For periodic boundary conditions 
in the direction of transfer, Z is given by the trace over this product and for free boundary 
conditions by a particular matrix element. Thermal expectation values can be calculated 
by inserting additional operators. In this way, the complete thermodynamical behaviour 
is contained in T. 
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Figure Bl. Sections of a square lattice (thick) which form (a) the row transfer matrix and (b) 
the corner transfer matrix. Open circles correspond to fixed spins. 



It is customary to write the transfer matrices as 

T = exp(-ff) (Bl) 

in analogy with time evolution operators. For the Ising model, T can be expressed in 
terms of Pauli matrices and the corresponding H may be regarded as the Hamiltonian for 
a certain spin one-half quantum chain. In the following this is described in some detail. 
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B.l. Row transfer matrix 
This quantity (Figure Bla) is appropriate for a homogeneous system or a layered one 
with couplings Ki(n) varying along the horizontal direction. It is given explicitely by [167] 



T 



exp 



exp 



(B2) 



where the K^(n) are dual couplings (tanhi^ = exp(— 2Ki)) and periodic boundary 
conditions are assumed horizontally, cr£ +1 = af . An open system can be obtained by 
setting one K 2 (n) equal to zero. 

The matrix T, or its symmetrized versions one uses normally, can in principle be 
diagonalized for arbitrary couplings [168, 78]. In this sense, the two-dimensional Ising 
model is called a completely integrable system. A considerable simplification occurs, 
however, in the anisotropic (Hamiltonian) limit of strong vertical and weak horizontal 
bonds [169, 170, 135]. Then K{(n) and K 2 (n) are both small and one can combine the 
exponentials to obtain H = 2K^H, where K{ is some reference value and 



H 



1 L 1 L 

~~ 9 Y kn(7 n ~ 9 Y A " a n<+1 



(S3) 



n=l 



This is the Hamiltonian of an Ising chain with couplings A n and transverse fields h n . The 
normalization is such that for h n = X n = 1 the excitations have velocity one. 

The operator H can be expressed in terms of Fermi operators c n , c+ via the Jordan- 
Wigner transformation [171, 167, 172] 



n-l 



r n = c+ exp 



1=1 



o~ n — 2c^c n 1 



(S4) 
(£5) 



This gives 

L 



L-l 



7< = ~Y hn \ C K Cn ~ \ I ~\ Yl A "( c n - c n)(c++i+c n+ i) + ^ A L (c£-c L )(c^ + ci)P (BQ) 



where V = (—^-) L Yln a n = ex P(' 7r E n c ii c ' 1 )- This operator with eigenvalues p = ±1 
commutes with 7i and distinguishes subspaces with even and odd fermion number. In each 
subspace H is diagonalized by a Bogoljubov transformation [80, 173]. In terms of new 
operators J it takes the form 

1' 



E q | Ci q ~OL q 



(B7) 
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The single-fermion eigenvalues e q follow from a L x L matrix equation which reads, in the 
notation of reference [80] 

(A-B)(A+B)* 9 = eJ*, (B8) 

with the matrix given by 

/ hj + Xl /11A1 -p h L \ L \ 



friAi hj + \( h 2 \ 2 



hL-2^L-2 h 2 L _ x + \\_ 2 hL-lXL-1 

\-ph L X L Hl-xXl-x h\ + \\_ x ) 



(B9) 



For a homogeneous system the equations are solved by Fourier transformation. The fermion 
states are running waves with momenta q given by q = 2nn/L (p = — 1) or q = (2n + 1)tt/L 
(p = +1) and integers n such that | q \ < tt. Putting h n = 1, A n = A, the eigenvalues are 

£q = (l + A 2 - 2Acosq) 1/2 . (BIO) 

This expression is the Hamiltonian limit of Onsager's general result [174] 

cosh uj q = cosh 2K{ cosh 2K 2 — sinh 2K* sinh 2K 2 cos q (511) 

where u q = 2K*e q . 

If one multiplies a large number of T's, their largest eigenvalues are most important. 
They are given by the smallest ones of H. The two lowest states of H are the ground states 
in the two subspaces p = ±1. If A > 1, they become degenerate for L — > oo and this leads 
to the appearance of long-range order in the system. 

If the system has free ends, no distinction occurs between p = ±1. The solutions 
are standing waves with g-values determined by the boundary conditions. For A > 1 one 
finds an additional solution with imaginary wave number, localized near one surface. The 
corresponding eigenvalue vanishes exponentially as L — > oo so that this state is again 
related to the long-range order. In particular, it determines the surface magnetization 
(cf. Section 5.2). 

At the critical (self-dual) point h = A = 1 the eigenvalues are 

e q = 2 | sin | | . (£12) 

The (/-values for the homogeneous case are given above, while for free ends q = (2n + 
l)n/(2L + 1) and < q < rc. The low-lying part of the spectrum has a linear energy- 
momentum relation, e q = q, and the towers of equidistant eigenvalues are in complete 
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agreement with the conformal predictions (Appendix A). The left- and right-moving 
particles in the homogeneous system correspond to the two Virasoro algebras present in 
this case. 

The layered systems treated in Sections 5-7 are special cases of the problem formulated 
in (B3) or (B8-9). 

For periodic boundary conditions the existence of two subspaces leads to considerable 
complications if an operator like a x connects the two. Thus a spin correlation function in 
the direction of the transfer has the form 

< ofo% >= \< m i "P I ffI >! 2 exp [-1 (Em(-p) ~ E (p))} (S13) 

m 

The non-universal behaviour of this function for the case of a line defect (Section 6) can 
be related to its particular structure involving both subspaces [119]. 

B.2. Corner transfer matrix 

This quantity is the partition function for a whole angular segment of a lattice. For 
the case shown in Figure Bib it is a 90° segment. The inner variables are supposed to be 
summed over. The partition function of the total system is then given by the trace over 
the product of four such 90° matrices. If the system is invariant under 90° rotations, all 
matrices are identical. Otherwise one has two different types. The importance of these 
matrices is connected with their use in calculating order parameters. For this the variables 
along the outer boundary are fixed as shown and the order at site is considered in the 
thermodynamic limit. That this is a practical and quite efficient procedure was noted 
by Baxter who introduced this type of transfer matrix and studied it for various solvable 
models [175, 176, 6, 177]. 

From Baxter's results it follows that the operator H related to the corner transfer 
matrix of a homogeneous Ising model has again the form (B3). In the Hamiltonian limit 
this can be seen directly. The coefficients h n , X n then follow from the geometry and are 
proportional to the number of vertical and horizontal bonds at distance n from the tip, 
respectively 

h n = 2n A n = A(2n + l) n = 0, 1, • • • (L - 1) (£14) 

In addition hi, = due to the fixed boundary spins. The resulting eigenvalue problem 
(B8-9) can be solved via generating functions and special polynomials [178, 179]. Due to 
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the boundary conditions, one eigenvalue is zero. The other low-lying eigenvalues in the 
ordered phase (A > 1) are, with q replaced by I 

ei = (2l-l)e Z = 1, 2, - • • (S15) 

where e = 7rA/2K(l/A) and K is the complete elliptic integral of the first kind. This result 
also follows by taking the Hamiltonian limit in the general formula for the eigenvalues 
uji = 2K*ei of the infinite system [6] 

where k = (sinh2i^i sinh2iir2) _1 and u, measuring the anisotropy, is defined via the elliptic 
function sn by sinh2.ftT2 = — isn(ht, k). The equidistance of the levels si, lo\ which directly 
reflects the geometry, is a very remarkable property. Such a level structure is also found 
for the corner transfer matrices of other integrable models. From (B16) it follows that near 
criticality the u)\ scale as l/ln(l/£).Thus the spectrum collapses at the critical point, but 
(B15) still holds for a large finite system. Then e = n/lnL in agreement with conformal 
predictions [35, 180, 181]. 

For radially inhomogeneous systems as in Section 8.2, the coefficients h n , A n are 
modified compared to (B14). The same holds if the shape is different from a wedge, 
as in Section 4.3. 

The magnetization < cr > in the centre of a system built from m segments and closed 
upon itself, can be written as 

< * >= t^fc < B17 > 

where Z + and Z- are the partition functions with ao parallel and antiparallel to the 
boundary spins, respectively. In terms of the corner transfer matrix T of one segment this 
becomes a quotient of traces 

The operator ctq(J% measures the relative orientation of the spin at and the outer spins. It 
can be expressed in terms of the Fermi operators which diagonalize T as exp(i7r ^ z otfai) 
where I > 1. Thus it distinguishes states with even or odd number of excited fermions. 
Inserting (Bl) and the diagonal form of H, the traces can be performed independently 
over the single-fermion states and one obtains the magnetization in the form of a product 

<a >=ntanh(^) (B19) 
l 
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This is a central result of the corner transfer matrix approach. 

One may note that the formula (B19) is identical to the one for the corresponding end 
magnetization of an inhomogeneous Ising chain with couplings mui/2. Onsager's result 
[182] for the spontaneous bulk magnetization < o~o >= (1 — A; 2 ) 1 / 8 follows directly from it 
by setting m = 4 and using (B16) together with elliptic function identities. As mentioned 
in Section 3.2, the situation is more complicated for systems having free edges [34]. 

B.3. Rescaling 




fa) (b) 



Figure B2. Rescaling of a system with different correlation lengths (a) into an isotropic one 
(b). The change of an angle is also indicated. 

The results of the Hamiltonian limit may also be used to discuss isotropic systems. 
The problem becomes effectively isotropic if the correlation lengths in the two directions 
are rescaled as shown in Figure B2. For this, one introduces lattice constants such 
that ai£i = a 2 ^2 [26, 31]. For the square lattice near criticality one has £2/^1 = 
cosh 2K 2 / cosh 2K t ^ 2K{ < 1. Thus a 2 = 1 gives a x = 2K{. The operator H for 
row transfer then corresponds to a unit step in the isotropic system. For the wedge shown 
in the figure, the original opening angle of 90° is reduced to 9 = AK{. Therefore the 
operator H for corner transfer corresponds to a step of twice the unit angle in the isotropic 
case. There is no change in the angle if the wedge is bounded by the principal axes. 
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Appendix C. Perturbation theory for extended defects 

A perturbation approach to the extended defect critical behaviour has been proposed 
by Bariev [141]. It provides a criterion for marginal behaviour and, when the system is 
marginal, allows a determination of local first-order corrections to the scaling dimensions of 
either the order parameter cr(r) or the energy density s(r). These corrections are obtained 
through a calculation of the two-point correlation functions to logarithmic accuracy, 
making use of the operator product expansion (OPE) [183, 184]. This is an extension 
to inhomogeneous sytems of techniques introduced by Kadanoff and Wegner [185] in 
statistical physics and Polyakov [186] in quantum field theory. In the following the method 
is presented for a temperature-like extended perturbation and the correction to the order- 
parameter correlation function is determined. 



C.l. General 

One considers, in the continuum limit, the following perturbed <i-dimensional system 



-PH = -PH + g 



J Z(r)e(r)dr (CI) 



where Hq is the unperturbed Hamiltonian and g the perturbation amplitude. The shape 
function Z(r), which gives the form of the inhomogeneity, is assumed to be homogeneous 
in r so that it may be generally written as 

Z(r) = M (C2) 

where u is a unit vector along r. The numerator contains the angular dependence of the 
perturbation which decays from the core of the defect with dimension d*. For a point 
source f(u) = 1 whereas for a line or a plane f(u) =| cos6> | _tiJ where 6 is the angle 
measured from a polar axis which is orthogonal to the core subspace. 

At the critical point of the unperturbed system the OPE is used, in multipoint 
correlation functions, to expand the product of two local operators at r\ and r 2 on the 
complete set of operators associated with the fixed point of Hq. The distance r±2 has to 
be smaller than the distances for other pairs appearing in the correlation function. Only 
the most singular term in the expansion is retained. The following reduction equations are 
needed: 

o-(n)e(r 2 ) ^ aa(r!)r^ 2 x e{r 1 )e{r 2 ) ~ be{r{)r^ (C3) 
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They are written for a translationally invariant system, x is the bulk scaling dimension of 
the energy density and averages are supposed to be subtracted. The powers in (C3) follow 
from scaling. 

In a semi-infinite system translational invariance is lost in the direction perpendicular 
to the surface. When r*i belongs to the surface, the structure constants a and b in (C3) 
only acquire an angular dependence. When r2 goes to the surface, i.e. when 9 ~ 7r/2 — rj^ 1 , 
the surface scaling dimension d (Section 2.3), has to replace the bulk one in (C3). Then, 
as for surface scaling functions [18], one expects quite generally a{9) ~ b(9) ~ (cos9) d ~ x . 

For the 2d Ising model with x = 1, one indeed obtains [187] 

a(9) = - cos 9 b(9) = - cos 6 sin 2 6 (C4) 

7T 7T 

The bulk values are a = l/2n and 6 = 0. The last one vanishes due to the dual symmetry 
(e — > — e) of the bulk Ising model in two dimensions [156]. 



C.2. Relevance-irrelevance criterion 
The order parameter correlation function has the following expansion in g 

G(R,g) = 9 —\ \ , « (T(0)a(R)e( ri ) ■ ■ -e(r n ) » Y[z(n)dri (C5) 

n=0 n ' ^ £ '' i=l 

where the double brackets denote the irreducible part of the multi-point correlation 
function. The integral extends over the subspace £' with dimension d' where the 
perturbation is nonvanishing. R is assumed to be orthogonal to the core subspace. 
The OPE allows an evaluation of the first-order correction in (C5) giving 

5G {1) ~ gG(R, 0)R d e"- x (C6) 

where x is the bulk scaling dimension of the energy density and, for a bulk perturbation, 
the effective dimension of the defect is defined as 

d e ff = max [d' - to, d*] (C7) 

The second alternative, d e ff = d*, is met when the angular integral in (C5) becomes 
divergent at 9 = ix/2 due to the singularity of Z(r) at the core of the defect. This is cured 
by introducing a cut-off at 9 = n/2 — r _1 which modifies the independence in (C6) [141]. 

When d e ff < x, the leading correction to the long-distance behaviour is small and 
the perturbation is irrelevant. In the opposite case, the correction factor in (C6) grows 
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ClclcCL lype 


LI 


LI 


U, 


deff 


narrow line defect (a; = 0) 


1 


1 


2 


1 


surface extended line defect 


1 


2 


2 


2-uj 


bulk extended line defect 


1 


2 


2 


max [2 — a;, 1] 


radial defect 





2 


2 


2-w 


surface field 





1 


2 


l-UJ 



Table CI. Values of the effective dimension of the defects, d e //, for the different systems 
studied in this review. 

as a power of R, indicating a relevant perturbation and the expansion becomes useless. 
Finally, when d e ff = x, the perturbation is marginal and logarithmic corrections leading 
to (7-dependent local exponents are obtained. As to the value of <i e //, Table CI lists it for 
the systems studied in this review. 

For slow decay (uj small), d e ff = d' — uj, the perturbation is long-ranged and marginal 
behaviour sets in when uj = d' —x. This is an extension to d! < d of the condition obtained 
in (5.2) via scaling arguments. For large enough uj, d e f / = d*, the perturbation is effectively 
short-ranged and one recovers the marginality criterion (6.5). Similar conclusions can be 
drawn for an order parameter perturbation as in Section 5.4. 

With a free surface, the angular integral is modified through the angular dependence 
of the structure constants in (C4). As a result, d* in (C7) has to be replaced by d* — d + x 
and marginal behaviour can no longer be induced by an effectively short-range defect since 
d* < d. 



C.3. Marginal behaviour 
In the following, one assumes a bulk temperature-like perturbation and a defect with 
d e ff = d' — uj = x > d* . The n-th order term in (C5) can be rewritten so that the 
integration extends over r\ < r 2 ■ ■ ■ < r n . The leading contribution to the multi-point 
function comes from pairs of points which are close to each other. Assuming that the 
shortest distance is between r*i and the origin and using the first reduction relation in 
(C3), the integral over r\ gives 

aa(0) ^ dnrf'- 1 -"-* J /(ui)du! (C8) 
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where we wrote dr\ = rf _1 dridn-i for the volume element. Since d e ff > d*, the angular 
integral is regular and contributes a geometrical factor Sd'd*(w) leading to 

» n 

SOW =g n aS d , d 4u) / lnr 2 « a(0)a(R)e(r 2 ) ■ ■ -e{r n ) > Y[Z(n)dri (C9) 

J£'(r 2 <r 3 ---<r n ) i=2 

The same process can be iterated n times giving 

SG {n) = ^ [<7a<WMhi#] n < a(0)a(R) > (CIO) 
n! 

It may be shown [141] that when other pairs like e(r 1 )e(r 2 ) are first contracted, the result is 
logarithmically smaller. Thus, to the leading logarithmic order, the perturbed correlation 
function is given by 

G(R, g) = G(R, 0)R 9aS *'<t* ( ^ (Cll) 

From it, one deduces the first-order local change of the bulk order parameter scaling 
dimension x, 

x^x-gaSd'd^^ + Oig 2 ) (C12) 
In the same way, for the local scaling dimension of the energy density, one obtains 

x x = x-gbS d , d 4u) + 0(g 2 ) (CT3) 

where x is the unperturbed bulk value [141]. 

These results are valid provided | g |<C 1 and | g \ In 2 R^> 1 and can be used in the scaling 
region. When the two alternative conditions for marginal behaviour are simultaneously 
fulfilled, i.e. when d' — uj = d* = x, the perturbation results are only valid in a crossover 
regime with | g \ lni? 1 and | g \ In 3 R ^> 1 [187]. Then, due to the onset of local order 
at the bulk critical point when g > 0, the local dimensions are singular in g at g = (see 
Section 7). 
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point 


source 


line 


source 


bulk order parameter 


1 

8 


-9 






surface order parameter 


1 

2 


7T f 


l 

2 




surface energy density 


2 - 


Us 


2 


-Ag 



Table C2. Local scaling dimensions for extended, temperature-like, marginal defects in the 2d 
Ising model, up to first order in the perturbation amplitude g. 

In the semi-infinite geometry, with d e ff = d! — oj = x, the calculation proceeds as 
in the bulk. The only change is introduced by the angular dependence of the structure 
constants [187]. Explicit results for marginal, bulk and surface extended perturbations in 
the two-dimensional Ising model are given in Table C2. 
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